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Abstract 

Single reference coupled-cluster (CC) methods for electronic excitation are based on a biorthog- 
onal representation (bCC) of the (shifted) Hamiltonian in terms of excited CC states, also referred 
to as correlated excited (CE) states, and an associated set of states biorthogonal to the CE states, 
the latter being essentially configuration interaction (CI) configurations. The bCC representation 
generates a non-hermitian secular matrix, the eigenvalues representing excitation energies, while 
the corresponding spectral intensities are to be derived from both the left and right eigenvectors. 
Using the perspective of the bCC representation, a systematic and comprehensive analysis of the 
excited-state CC methods is given, extending and generalizing previous such studies. Here, the 
essential topics are the truncation error characteristics and the separability properties, the latter 
being crucial for designing size-consistent approximation schemes. Based on the general order 
relations for the bCC secular matrix and the (left and right) eigenvector matrices, formulas for 
the perturbation-theoretical (PT) order of the truncation errors (TEO) are derived for energies, 
transition moments, and property matrix elements of arbitrary excitation classes and truncation 
levels. In the analysis of the separability properties of the transition moments, the decisive role 
of the so-called dual ground state is revealed. Due to the use of CE states the bCC approach can 
be compared to so-called intermediate state representation (ISR) methods based exclusively on 
suitably orthonormalized CE states. As the present analysis shows, the bCC approach has decisive 
advantages over the conventional CI treatment, but also distinctly weaker TEO and separability 
properties in comparison with a full (and hermitian) ISR method. 

PACS numbers: 
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I. INTRODUCTION 



The extension of the coupled-cluster (CC) approach PQI21E], originally devised for ground- 
states, to the treatment of electronic excitation has afforded the emergence of a variety 
of highly successful computational methods, excelling at the potential for both numerical 
efficiency and accuracy. The excited-state CC methodology comprises three major devel- 
opments referred to as coupled-cluster linear response (CCLR) theory IH El [6l [71 E], the 
equation-of-motion coupled-cluster (EOM-CC) approach [9l [TOl [TTl [121 [131 E], and the 
symmetry-adapted cluster configuration interaction (SAC-CI) [ISl [ISl [IZj ■ While these de- 
velopments vary in the derivation of the CC equations, the resulting computational schemes 
are largely equivalent. A notable difference, however, is the treatment of transition moments 
and excited-state properties, where, in contrast to the CC-EOM and SAC-CI schemes, CCLR 
theory leads to size-consistent expressions. 

The basic feature of the CC methods is the use of so-called correlated excited (CE) states 
as basis states in the expansion of the exact excited states. These CE states are obtained 
by applying physical excitation operators associated with single (S), double (D), triple (T), 
. . . electron excitations to the exact (correlated) CC ground-state, |\E'j) = Cj\'^q'^), rather 
than to the Hartree-Fock (HP) ground state, = Cj\^o), establishing in the latter case 
the familiar configuration interaction (CI) basis states or CI configurations. Obviously, the 
CE states represent intermediates of sorts, positioned in-between the simple CI configura- 
tions and the exact final states. Accordingly, methods based on the use of CE states have 
been referred to as intermediate state representations (ISR) [18J, and the excited-state CC 
approach is closely related to the family of ISR methods. However, the CE states are not 
orthonormal, and this problem is dealt with by introducing in addition to the CE states the 
associated set of biorthogonal states. Using the biorthogonal states on the left side and the 
CE states on the right-hand side, one obtains a mixed (biorthogonal) representation of the 
(shifted) Hamiltonian giving rise to a non-hermitian secular matrix. The two sets of states 
differ distinctly in their intrinsic quality. In fact, the biorthogonal states can be identified 
essentially as CI configurations. This means that the bCC (biorthogonal CC) representation 
represents a hybrid approach, combining the CI and ISR concepts in equal measure. 

What are the merits of an ISR approach as compared to the conceptually so much sim- 
pler CI treatment? The answer is that the ISR methods are not (or much less) affected by 



3 



two basic deficiencies of tlie CI approximation scliemes, namely tlie lack of size-consistency 
and the need for relatively large explicit configuration spaces. The size-consistency error 
inherent to limited CI treatments (as opposed to full(F) CI expansions) stems from the non- 
separable structure of the CI secular equations. More precisely, for a system composed of 
non-interacting fragments there is no a priori decoupling of the CI secular matrix into cor- 
responding fragment blocks. The inevitable and, moreover, uncontrollable size-consistency 
error associated with truncated CI expansions grows with the system size, rendering the 
results for extended systems useless. This is why CI cannot rank as a genuine many-body 
method. Secondly, large CI expansions (configuration spaces) are needed to suppress the 
quite unfavorable truncation error, that is, the error due to discarding higher excitation 
classes in the CI expansions. For example, a CI expansion extending through single and 
double excitations (CISD) induces a truncation error of second-order of perturbation theory 
(PT) in the excitation energies of singly excited states. In the CC methods, by contrast, the 
corresponding (SD) truncation error is of third order. Yet more advantagous is the situation 
in full (hermitian) ISR methods, where the SD truncation error is already of fourth order. 

The purpose of this paper is to review the separability properties and truncation error 
characteristics of the excited-state CC schemes, aiming here at a more systematic and more 
comprehensive analysis than available so far. Previous studies of this kind have been pre- 
sented by J0rgensen and collaborators within the context of CCLR theory, addressing the 
size-consistency of the CC excitation energies [S] and transition moments |T9J, and analyzing 
the truncation errors in the CC excitation energies [201 EI]- Studies devoted to the size- 
consistency of the excited state CC equations have also been presented by Mukhopadhyay 
et al. |22] and Stanton [23] . From a different perspective, the so-called order relations of the 
bCC secular matrix and the ensuing truncation errors in the CC excitation energies have 
been discussed by the present authors [18] and by Trofimov et al. [23]. However, the latter 
studies, devised essentially to enable comparison to the ADC (algebraic-diagrammatic con- 
struction) propagator methods [2ll l25l [26] , were somewhat limited and suffered, moreover, 
from certain misconceptions concerning the CC transition moments. 

The analysis given in this paper of the excited-states CC approach will be formulated 
entirely within the framework of the bCC concept, that is, a non-hermitian secular problem 
associated with a dual representation of the (shifted) Hamiltonian in terms of two biorthog- 
onal sets of basis states. While in the EOM-CC development, the bCC representation is 
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introduced more or less in an ad hoc manner, depicting the bCC secular matrix essentially 
as a Cl-type representation of an effective (similarity transformed) Hamiltonian, the CCLR 
approach is based on response theory for time-dependent CC ground state expectation val- 
ues. Clearly, the CCLR derivation is highly original and instructive, and, in fact, transcends 
the simple bCC formulation in the case of the transition moments and various reponse 
properties. However, the essence of the emerging computational scheme can much easier 
be presented and understood using the bCC concept. Thus, it should be permissible and 
even advisable to abandon the original notations associated with linear response theory and 
rather resort to a stringent wave function formulation adapted to the bCC concept. Not only 
will this make the excited state CC methods more amenable to readers not familiar with 
the rather demanding time-dependent CC response theory, it will also allow us to embed 
the CC approach quite naturally in the broader context of ISR methods. 

An outhne of the paper is as follows. In Sec. II, we briefly review the CI method with 
regard to the truncation error and separability properties. Sec. Ill presents the basics of 
the CC approach to electronic excitation. This is followed in Sec. IV by an analysis of 
the properties of the bCC representation and the resulting excitation energies, transition 
moments, and excited state properties. Sec. V contrasts the bCC representation with a 
full (hermitian) ISR approach. A summary and some conclusions are given in the final 
Sec. VI. Some important supporting material is presented in a tripartite Appendix. In 
App. A the proof of the bCC order relations is reviewed. In App. B we derive the order 
relations of both CI and bCC eigenvector matrices, which, in turn, allows us to generate 
general truncation error formulas. The CCLR forms of the right transition moments and 
excited-state properties are reviewed in App. C. 
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II. A LOOK AT THE CI METHOD 



The justification of tlie intrinsically more complicated intermediate state representations 
derives from basic shortcomings of the standard configuration-interaction (CI) approach 
with regard to the truncation error of the CI expansions and the size-consistency of the 
results. To provide for a general background we begin with a brief recapitulation of the CI 
problems. 

For the electronically excited states of an atom or molecule the Schrodinger equation 
may be written in the form 

(^-Eo)|^n) =^n|^n) (l) 

where H is the Hamiltonian of the system under consideration, and Eq and ujn = En — Eq 
denote the ground state energy and excitation energy, respectively. In the CI treatment, 
being the standard quantum chemical method, the excited states are expanded according to 

|$„) = ^Xj„|<I>j) (2) 
J 

linear combination of CI states 

\^j)=Cj\^,) (3) 

generated by applying "physical" excitation operators Cj to the Hartree-Fock (HF) ground 
state, |$o)- Using the notation of second quantization, the excitation operator manifold may 
be expressed as follows: 

{Cj} = {c\ck] clclckQ, a<b,k <l;...} (4) 

Here cj,(cp) denote creation (annihilition) operators associated with HF orbitals Follow- 
ing a widely adopted convention, the subscripts a,b,c, . . . and k, . . . denote unoccupied 
(virtual) and occupied orbitals, respectively, while the indices p,q,r, . . . will be used in 
the general case. The capital indices I,J,... are used as an abbreviation for strings of 
one-particle indices, e.g., / = (abkl). 

The excitation operators in ^ can be divided into classes of p-h (single), 2p-2h (double), 
. . . excitations. For brevity, these classes will be numbered consecutively, that is, the class 
of fip-fih excitations is referred to as class fi. The class of a particular excitation J will be 
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denoted by [J]. For example, [J] = 2 means J is a double excitation. The HF ground state 
|$o); being part of the CI expansions ([2]), constitutes a zeroth class {fi = 0). 

The excitation energies and expansion coefficients are obtained as the eigenvalues and 
eigenvector components, respectively, of the CI eigenvalue problem, reading in a compact 
matrix notation 

Hx = xn, X^X = 1 (5) 
Here H is the (subtracted) CI secular matrix, 

Hjj={^j\H-Eo\^j) (6) 

fl denotes the diagonal matrix of excitation energies Un, and X is the matrix of (column) 
eigenvectors X„. Note that the subtraction of the ground-state energy Eq in the diagonal of 
the CI secular matrix is a mere convention here, introduced for formal analogy to the bCC 
representation considered in Sec. III. 

Approximate CI treatments are obtained by limited CI expansions as opposed to full 
(FCI) expansions. In the following we will be concerned with systematic truncations of the 
CI expansions, that is, expansions being complete through a given excitation class fi. These 
systematic truncation schemes can be examined with respect to the perturbation-theoretical 
order of the induced error in the CI results. For this purpose one has to inspect the "order 
structure" of the CI secular matrix (Fig. 1), that is, a perturbation-theoretical classification 
of the subblocks Hfj,^, associated with a partitioning of H with respect to the excitation 
classes. Fig. 1 shows the characteristic CI structure, where each excitation class is coupled 
linearly in the Coulomb integrals (ffist order) to the next and next but one excitation class. 
Owing to the (zeroth order) orbital energy contributions of the diagonal matrix elements, 
the diagonal blocks are indicated by zeros. 

The order structure of the CI secular matrix gives rise to characteristic truncation errors 
in the excitation energies. As the most important case, let us consider singly excited states, 
that is, states deriving perturbation-theoretically from p-h configurations (CI basis states). 
Due to the linear (ffist-order) coupling to triple {3p-3h) excitations (see Fig. 1), there is a 
second-order energy contribution to the single excitation energies arising from the admixture 
of triple excitations. This means that a second-order truncation error arises in the single- 
excitation energies if the triple excitations are not taken into account. A stringent derivation 
of the general CI truncation errors is given in App. B.l. Specifically, the truncation error 
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orders (TEO) for single excitation energies are given by the following formula: 

{u, a even 

(7) 
+ 1, fjL odd 

Here n denotes the highest excitation class included in the CI expansion manifold. In Table 1 
the resulting TEOs for = 1, . . . , 6 are listed. 
Besides the energies, the transition moments 

7; = (*„|Z)|*o) (8) 

are of interest, required to compute spectral intensities. Here D denotes a (one-particle) 
transition operator, e.g., a dipole operator component. To evaluate the truncation error of 
the transition moments one has to analyze the CI expression 

Tn = KiDX, (9) 

with respect to the order relations of the CI eigenvectors X„ and Xq, respectively, and the 
CI representation of £>, 

Dij = ($7|D|$j) (10) 

The order relations of the CI eigenvectors for singly excited states and for the ground state, 
shown in Fig. 1, are part of the general order structure of the CI eigenvector matrix X, 
derived in App. B.l. Together with the trivial order structure of D (zeroth-order diagonal 
and off-diagonal blocks, other matrix elements vanishing) one can deduce the following 
simple TEO expression for the transition moments of singly excited states: 

Ote{i^) = /i (11) 

For odd values of the TEOs of the transition moments are smaller by 1 as compared to 
the excitation energies. 

The CI truncation errors arc relatively large, which, in turn, implies that large CI ex- 
pansions are required to meet specific accuracy levels. For example, in order to treat singly 
excited states consistently through second order of perturbation theory (PT), the CI configu- 
ration space must comprise the triple excitations {/i — 3). By contrast, in the ISR methods a 
much smaller explicit configuration space, consisting of single and double excitations, affords 
the corresponding level of accuracy. 
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To analyze the size-consistency properties of a method, one usually resorts to the separate 
fragment model, that is, a system 5* consisting of two strictly non-interacting fragments, A 
and B. A method for treating electronic excitation is size-consistent (here, more specifically, 
size-intensive) if for local excitations, say on fragment A, the computed excitation energies 
and transition moments do not depend on whether the method is applied to the fragment 
or the composite system. It is well known that truncated CI treatments do not fulfill 
this property. There is an uncontrollable size-consistency error in the treatment of the 
total system, corrupting not only the results for the separate fragment model, but, more 
generally, any truncated CI treatment of larger molecules. This can be nicely demonstrated 
in the exactly solvable model of a chain of non- interacting two-electron two-orbital (2£'-20) 
systems, such ls^2s° He atoms or minimal basis H2 molecules (see, for example, Meunier 
and Levy [27]). Because of this deficiency the CI method does not qualify for a genuine 
many-body method. 

As a preparation for the analysis of the bCC schemes let us briefly inspect the case of CI in 
some more detail. Obviously, the Hamiltonian for S decomposes into the sum of the fragment 
Hamiltonians, H = Ha + Hb- Moreover, the one-particle states (HF orbitals) of S can be 
classified as belonging either to fragment Aoi B (local on A or B, respectively). Accordingly, 
the CI states can be partitioned into three different sets, that is, local excitations Ia on 
fragment A, local excitations Ib on fragment B, and mixed (or non-local) excitations Iab 
involving both fragment A and B. In the latter class we may disregard any excitations that 
do not conserve the local electron number, e.g., A^B^ charge-transfer excitations. In the 
non-interacting fragment model, such charge transfer excitations are strictly decoupled from 
the fragment-charge conserving excitations, which are of interest here. 

The CI configurations 1$/) can be written as products of fragment states, e.g.. 



It should be noted that the neglect of full antisymmetrization of these product states is 
irrelevant because in the non-interacting fragment model matrix elements are not affected 
by inter-fragment antisymmetrization. Fig. 2 shows the partitioning of the CI secular ma- 
trix with respect to the three different types of configurations, that is, local excitations on 




$0) 



*o)lO 



(12) 
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fragment A, local excitations on fragment B, and non-local excitations, respectively. While 
the A and B states are strictly decoupled, Hab = 0, there is a coupling between the local 
and mixed excitations, e.g., 

Hj,,j,j,=6j,jAK\Hb\^%) (13) 



as can easily be derived using the CI state factorization according to Eq. (12). An explicit 
example is the coupling matrix element 

Hak,bc'd'i'j'l = ^ab^klVc'd'li'j'] (14) 

for a single excitation = {ak) associated with fragment A and a non-local triple excitation 
Jab = {bl,c'd'i'j'). Here the unprimed (primed) indices denote fragment A (fragment B) 
one-particle states; Vpg[rs] denotes the anti-symmetrized Coulomb integral. 

Given the structure as shown in Fig. 2, the CI secular matrix is said to be non-separable 
because there is no a priori decoupling of local excitations (say on fragment A) from non- 
local (or mixed) excitations. The CI treatment of the composite system 5* aims in an 
inextricable way at an optimal description of both fragments, that is, the excited state of 
fragment A and the ground state of fragment B. In the exact (full) CI result, say for the 
energy En = + Eq of the locally excited system, the ground state energy Eq of the 
unaffected fragment B would cancel exactly upon subtraction of the exact ground-state 
energy Eq = Eq + Eq of S, so that, of course, the full CI excitation energy. En — Eq = E^ — 
Eq, is size-consistent. At the level of a limited CI expansion, however, neither the excited 
state energy nor the ground-state energy are simply the sums of the fragment energies. 
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III. BIORTHOGONAL CC REPRESENTATION 



The bCC formulation (see, for example, Helgaker et al. |2H]) is based on a mixed repre- 
sentation of the (subtracted) Hamiltonian, H — Eq, in terms of two sets of left and right 
expansion manifolds: i) the CC states 

|vl>o>=C'j|vl>-) = C'^e^|$o) (15) 

on the right-hand side, and ii) the associated biorthogonal states 

($,| = {%\C]e-^ (16) 



on the left-hand side. Here Ci denote the physical excitation operators as specified in Eq. (|4 
The familiar ground state CC parametrization (and normalization) 

\n)=e^\^o) (17) 

is used, where |$o) denotes the HF ground state, and 

f = J2tiCi (18) 

is the cluster operator with the amplitudes tj determined by the ground-state CC equations. 
The T operator, comprising physical excitation operators only, commutes with any (physical) 



Cj operator, so that the CC states of Eq. (15) can likewise be written as 



|vI/0> = e^C'j|$o) (19) 

The biorthonormality of the two sets of states, 

($,|^°> = ($,| e-^e^ = Sjj (20) 

is an obvious consequence of the orthonormalization of the CI configurations = Cj |$o)- 
The bCC representation of H — Eq gives rise to a non-hermitian secular matrix M with 
the elements 

Mij = {^i\H-Eo\^^j) 

= {<^o\Cje"^[H,Cj]e'^\%). (21) 
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In the latter form the (CC) ground-state energy Eq no longer appears explicitly. Let us note 
that the bCC secular matrix can also be written as a CI representation 

Mjj = {<^i\H - Eo\<l>j) 

of the similarity transformed Hamiltonian, H = e~^He^. However, this form is less trans- 



parent then the bCC representation (21) and, thus, less useful for the analysis intended 
here. 

The (vertical) electronic excitation energies, Un = En — Eq, can be identified as the 
eigenvalues of the CC secular matrix M. Because M is not hermitian one has to deal with 
right and left eigenvalue problems, 

MX = Xri (22) 
Y^M = riY^ (23) 

where fl is the diagonal matrix of eigenvalues Un, and X and Y denote the matrices of the 
right and left eigenvectors, respectively. To obtain a definite normalization the two sets of 
secular equations have to combined according to 

Y^MX = ri, Y"fX = l (24) 

so that the resulting right and left eigenvectors form two mutually biorthonormal sets. As 
a consequence, the corresponding right and left excited states, 

1^-) = 5^X,„|vI/?> (25) 

(^^i^l = (26) 

are biorthonormal too, (\l/m = ^mn- 

In general, the right excited states I^I/Jf) are not yet eigenstates of H, because the underly- 
ing expansion manifold of Eq. (15) is incomplete as long as the CC ground state I^I/q'^) is 
not taken into account. Using the extended CC expansion manifold {R} = {I^E'o^) , |^/)} on 
the right-hand side, and, likewise, the extended biorthogonal manifold {L} = {($o| ? ("^/Il 
on the left side, one arrives at the full bCC representation of H — Eq associated with the 
extended secular matrix 



M' = I - I . (27) 
M 
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Here is a transposed (row) vector with the elements 

VI = ($01^1*?) (28) 

that is, the coupling matrix elements between the HF ground state and the excited CC 
states. Let us note that in the usage of the EOM-CC approach M' is denoted by H; in the 
CCLR context, on the other hand, the (inner) bCC secular matrix M (Eq. [2l| ) is referred 
to as the CC Jacobian A and the coupling vector is denoted by r]. 

The two expansion manifolds used in the bCC representation are of quite different quality. 
The "correlated excited states" (CES) of the set {R} are superior to the biorthogonal {L} 
states, if more complex. Obviously, a CC state of class [/] can be written according to 

|M/?>=e^|$,) = |$,)+ ""kI^k) (29) 

K, [K]>{I] 

as a linear combination of 1$/) and CI configurations of higher excitation classes, [K] > [I], 
extending through A^-tuple excitations. By contrast, the CI expansion of a biorthogonal 
{L} set state reads 

($,| = (<l>,|e-'^=($,|+ Yl (30) 

K, [K]<[I] 

that is, a linear combination of ($/| and lower class CI excitations, [K] < [I], including the 
zeroth class, [K] = 0. This follows from the observation that z^'* = ($/| e^'^ \ ^k) vanishes 
for [K] > [I] (and = ^ik for [K] = [/]). 

As is easily seen, the linear space spanned by the biorthogonal states through a given 
excitation class /i is identical with the corresponding space of CI configurations: 

span{{<!>i\ e"^, [/] = 0, 1, . . . , /x} = span{($/| , [J] = 0, 1, . . . , /i} (31) 

This means that (truncated) expansions in terms of the biorthogonal {{L} set) states are 
essentially of Cl-type. Let us note that the SAC-CI equations are obtained by using the CI 
expansion manifold for the left eigenstates rather than the {L} states [THl HSl [T7] . 

As will be discussed below, the use of a Cl-type expansion manifold on the left-hand side 
of the bCC representation deteriorates the overall order relations and separability properties 
to a certain extent. One may wonder then why one could not simply use the CC states as 
the common expansion manifold for both sides of the secular matrix. However, in such an 
approach, referred to as variational or unitary CC version (see Kutzelnigg [29j and Szalay 
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et al. [30]), there is a major problem associated with evaluating the secular (and overlap) 
matrix elements: being of the form ($/| e^'' He^ there is no obvious truncation of higher 
excitation contributions (below A^-tuple excitation level). 

Let us now briefly inspect the eigenpair manifold of the extended bCC secular matrix 
M'. Obviously, there is one more eigenvalue, ujq = 0, corresponding to the ground state, 
while the excited state eigenvalues u;„, n > 0, of the subblock M are also eigenvalues of M'. 
The corresponding extended left and right eigenvectors can easily be determined. Let us 
first consider the ground state solutions. Here the right eigenvector is trivial. 



^0 = (32) 




which is consistent with the fact that I^E^q^) is the exact ground-state. Less obvious is the 
left ground-state eigenvector 

= (l, Yl) (33) 
where the row vector can explicitly be obtained from M and v according to 

yI = -v'M-^ (34) 

The corresponding representation of the ground state, 

(^o| = ($o| + $^>^;o($/| (35) 

in terms of the biorthogonal bCC states is referred to as the "dual" ground state [7] . In the 
CCLR nomenclature, the dual ground state is denoted by (A|. As will be discussed below, 
the dual ground-state is a non-separable Cl-type representation of the ground state, leading 
to undesired features in the bCC transition moments. 

The left excited state eigenvectors of M' are obtained as obvious extensions according to 

ri = (0, n) (36) 

from the left eigenvectors of the M subblock. This means that the excited states {'^n'^\ of 



Eq. (26) are proper eigenstates of H. In particular, they are orthogonal to the exact CC 
ground-state, 

(vI/«|vI/-)=0 (37) 
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which follows from ($/|\&q'^) = 0. 

In general, that is, if not forbidden by symmetry, the right extended eigenvectors acquire 
non- vanishing zeroth components Xn = X'^q, and the extended eigenvectors take on the form 

where X„ is the (n-th) eigenvector of the M subblock, and the zeroth component x„ is given 
by 

a;„ = a;-VX„ (39) 
Since MX„ = uJnX_n the following relations hold: 

Xn = v'yi~^X^ = -Y}^X^ (40) 



Here Eq. (34) has been used to arrive at the last expression. 



As a consequence, the right expansion of an excited eigenstate takes on the form 

|v^W>=xJvI;-) + |vl/-) (41) 



where \^'^n) is given by Eq. (25). Let us note that (<l>o|^Jf) = since ($o|^/) = 0, so that 
the relation 

Xn = (^ol^i''^) (42) 

can be established. The excited eigenstates Y^n'^) are manifestly orthogonal to the dual 
ground state: 

(Wo|vl/M> = x„ + i:Sx„ = (43) 



where the relations = X/„ and Eq. (40) have been used. 

For spectral intensities the squared moduli |T„p of the transition moments ([s]) are re- 
quired, involving normalized ground and excited states. In the bCC representation a prop- 
erly normalized expression for |T„p is obtained according to im |19] 



|T„|2 = (^o| D {¥^\b 1*^^) (44) 

using both the left and right transition moments, 

= (45) 
tM = (^ol^l^i^)) (46) 
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Individually, the left and right transition moments have no significance because the respective 
ground and excited states are not normalized. However, the biorthonormality relations 



ensure the combined normalization in the product i\AAh. In 



contrast to the left transition moments, the ordinary bCC form (46) for the right transition 
moments is not separable, so that the results obtained at truncated bCC levels are not 
size-intensive [19]. Within the CCLR framework this shortcoming is avoided, as here a 
separable, if more elaborate expression is employed for the right transition moments [71 IT9] . 
The different treatment of the spectral intensities is a distinguishing feature of the otherwise 
equivalent CCLR and EOM-CC methods. In Sec. IV.D and App. C the CCLR expression 
for the right transition moments will briefly be reviewed. 

As will be discussed below, the problem in the right transition moments does not result 
from the right eigenstates but rather from the use of the dual ground state (^E'ol- As can be 



concluded from Eqs. ( [30|34[35[ ), the dual ground state is a Cl-type expansion of the form 



K>0 



where the expansion coefficients depend (via Eqs. 34p5 ) on the t-amplitudes of the CC 



ground state. For the CC ground-state energy the following relation holds: 

E^'^ = {%\H\%'^) (48) 

It should be noted that this expression applies not only to the exact CC ground state (where 
Eq'^ = Eq), but also to CC approximations based on truncations of the expansion manifolds. 



such as in CCSD. This can be seen by writing the rhs of Eq. (48) more explicitly as 



(vl/ol H \%') = ($01 ^ m) + J2 ^/o H (49) 

The first term on the rhs is the CC energy equation, while the second (summation) term van- 
ishes in compliance with the CC amplitude equations, ($/| H \'^'^) = ($/| e~^He^ |$o) = 0. 

What is the relation of the dual ground state to the CI ground state in the case of 
truncated expansions? Obviously, the energy expectation value of the dual ground state will 
always be greater than (or equal to) the corresponding CI energy. 



> E^' (50) 
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because the dual state expansion coefficients are non-variational. This can be nicely demon- 
strated in the exactly solvable model of 2 (or more) non-interacting 2E-20 systems (He 
atoms, or H2 molecules). 

It should be noted that the right and left eigenvalue problems ( 22p3 ) for the secular 
matrix M (Eq. 21) follow from a variational principle, = 0, under the constraint 

^$|^)> = 1. Independent variations and on the left and right side of the energy and 
overlap matrix elements lead directly to the right and left eigenvalue equations, respectively. 
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IV. ANALYSIS OF BCC EXCITATION ENERGIES AND INTENSITIES 



A. Order relations and separability of the bCC secular matrix 

Fig. 3 shows the order structure of the bCC secular matrix M. For the partitioning 
according to excitation classes, /i = 1, 2, . . . , the lowest (non-vanishing) PT orders are given 
here in the respective M^j, subblocks. In the upper right (UR) triangular part we recover the 
characteristic CI structure of Fig. 1. This outcome can readily be understood by inspecting 
the general expression for the secular matrix elements, 

($,| H |vl/0) = ($,| H 1$,) + J2 H 'k4'^ {^k\ H 1$^) (51) 

[K]<[I] [L]>[J] 

obtained by using the expansions ( 29p0 ) for the CC and biorthogonal states, respectively. 



For the UR matrix elements with [/] < [J], the sums on the rhs of Eq. (51 ) do not contribute, 
because the excitation classes of the double summation indices K and L differ at least by 
a triple excitation, [L] — [K] > 3, so that the Hamiltonian matrix elements {^k\ H 
vanish. This means that the bCC and CI secular matrix elements are identical, Mjj = Hjj, 
for [/] < [J]. For the diagonal blocks M^^, of course, the lowest non- vanishing PT order is 
zero because the perturbation expansions of the diagonal matrix elements Mjj begin with 
zeroth-order (HF) excitation energies. 

By contrast, the lower left (LL) triangular part, [/] > [J] gives rise to the remarkable 
"canonical" order relations, reading 

Mjjr^Oi[I]-[J]), [I]>[J] (52) 

This means that the lowest non-vanishing contribution in the PT expansion of the matrix 
element M/j, [/] > [J] is of the order [I] — [J]. Likewise we will use the notation 0[Mjj] = 
[I] — [J] . These order relations were first specified by Christiansen et al. [31] , quite explicitly, 
for the lowest 5 excitation classes (singles through pentuples) and later by Hald et al. [21j 
for general levels of excitation. A first general proof of the bCC order relations was given in 
Ref [TS]. A brief recapitulation of this proof is given in App. A. 

It should be noted that these canonical order relations are highly non-trivial indeed. Let 
us consider, for example, the niatrix elements, being of the order 2. This means that 



the apparent first-order contribution arising from the leading (CI) term on the rhs of Eq. (51 ) 
is exactly cancelled by other first-order contributions from the summation part. 
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The bCC order structure gives rise to specific truncation errors of tlie excitation energies, 
whicli can be analyzed by inspecting formal perturbation-theoretical (PT) interaction paths 
on the subblocks of the bCC secular matrix. As an example, let us consider the following 
path: 

Mn(0) ^ Mi3(l) ^ M33(0) ^ M3i(2) ^ Mn(0) 

This path allows one to specify formally a PT contribution in the excitation energy of a 
single excitation (class 1) arising from the admixture of triple excitations (class 3). The 
PT order of this path is 3, and there is no lower-order path involving class 3. This means 
that the PT contributions to single excitation energies arising from triple excitations are 
of the order 3. Stated differently, the truncation error for single excitation energies due to 
the neglect of triple excitations in the explicit bCC expansion manifold is of third order, as 
compared to a second-order truncation error in the CI treatment. 

A more general and stringent derivation of the bCC truncation errors is to be based on 
the order relations of the bCC eigenvector matrices, as is discussed in App. B. In short, 
the order relations in the bCC secular matrix induce corresponding order relations for the 
left and right eigenvector matrices Y and X, respectively. As shown in Fig. 10, the bCC 
eigenvector matrices combine canonical and Cl-type behaviour: the LL part of X and the 
UR part of Y exhibit canonical order relations, whereas Cl-type order relations (Fig. 9) 
apply to the UR part of X and the LL part of Y. The eigenvector order relations, in turn, 
allow one to analyze the truncation errors in the excitation energies, transition moments, 
and excited state properties. For the singly excited states, the truncation error orders (TEO) 



in the bCC excitation energies are given by the following formula (deriving from Eq. B.32) 



l/i, u even 
OTEif^) = r (53) 
+ /i odd 

Here n denotes the highest excitation class in the expansion manifold. In Table 1, the TEOs 
in the excitation energies of singly excited states are listed for the first 6 truncation levels. 

To discuss the separability properties of the bCC schemes [HI [191 EH] (see also Refs. [22], 
23] ) we revisit the separate fragment model, S = A + B, considered in Sec. 2. Like the 
Hamiltonian, H = Ha + Hb, also the CC operator can be written as the sum of the fragment 
operators, T = Ta + Tb- As a consequence, both the CC states |\E'5) the biorthogonal 
states ($/| used in the right and left expansion manifolds, respectively, can be written as 
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products of fragment A and B states: 

\K) = \K)\^o), \KJ = \K)\K) 

{^i. I = (K I ("^o"^ I ' i^iAis I = (K I (K I (54) 

Here the notation of the fragment states has been somewhat simphfied by omitting the 
superscripts and cc: for example, I^E';^) = "^ff^ l^^o) = l^^'*^^)- -^ig. 4 shows 
the partitioning of the bCC secular matrix with respect to the three different types of 
configurations, that is, local excitations on fragment A, local excitations on fragment B, 
and non-local (or mixed) excitations. Let us note once again that in the latter set we can 
disregard any charge-transfer type configurations. 

The separability properties of the bCC secular matrix, as shown in Fig. 4, can be readily 



derived by using the factorization (Eqs. 54) of the bCC basis functions. As an explicit 



example let us derive that the block Mab,a vanishes: 

= (KlHAl^jJiKl^o) + (Kl'^jJiKlHBl^o) = (55) 



Here the first term on the rhs vanishes because of the orthogonality relation ($j^|\E'^) = 0. 
The second term vanishes, 

{K\Hb\^o) = E^{K\'^o) = (56) 

because I^E'^) is an eigenfunction of Hb, and ($/^| and {"^q) are orthogonal. In fact, this 
result is not predicated on the exact fragment B ground state, but applies also to any (sys- 
tematic) CC approximation, since {^j^\Hb\'^q) = is satisfied as CC amplitudes equation 
for fragment B. 

In a similar way one may establish that the bCC secular matrix for fragment A is identical 
to the {A, y4)-block of the composite system secular matrix, 

Maa = (57) 

In deriving this result, one utilizes the CC equation for the ground-state of fragment B, 
{^q\Hb — EqI'^q) = 0. The matrix elements in the non-vanishing {A,AB) coupling block 
read 

MiMs = KjA^t (58) 
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where f = Finally, the non-local diagonal block matrix elements are given 

by 

Mui,,j,j, = 6j,j,Mf^j^ + 6i,j,mIj^ (59) 

It should be emphasized once again that these results do not presuppose exact CC ground 
states but apply as well to the (systematic) CC approximations. 

The separability structure of M reflects once more the different quality of the left and 
right bCC expansion manifolds. While the LL triangular is separable, the UR triangular 
displays the non-separable Cl-type structure of Fig. 2. What are the consequences for 
excitation energies and eigenvectors? Notwithstanding the apparently non-separable block 
structure of the bCC excitation energies are obtained in a separable way [8] . It is readily 
seen that the characteristic polynomial for the fragment secular matrix is a factor in the 
full characteristic polynomial associated with M. Accordingly, the eigenvalues (excitation 
energies) of fragment A are a subset of the eigenvalues of the full secular matrix M. This 
means that the energies of local excitations are separable quantities: the bCC results do not 
depend on whether the method is applied to the fragment or the composite system. 

For the eigenvectors one will expect different separability properties in the left and right 
manifolds. In fact, the right eigenvectors are separable, as shown in Fig. 4. For a local 
excitation n, say on fragment A, the only non-vanishing components are fragment-A com- 
ponents ^Laui ^'^d since TVf'^X^^ = ujn 2LAn is readily established that the fragment- A part 
of 2Ln is equal to the corresponding fragment-A eigenvector, 2LAn — 2Ln- More explicitly, 
the separability properties of right eigenvectors may be written as 

^lAn = XfAn-i ^Ibu = ^Iabu = 0, (60) 

assuming here a fragment- A excitation, n = ua- 

For the left eigenvectors y„ the fragment-A eigenvector is recovered by the A part of 
the full eigenvector, = Y_^. The local fragment-i? components vanish, Y/^„ = 0, but 
there are non-vanishing non-local components, Yj^^^n 7^ 0. These non-local components are 
related to the local ones according to 

Y-Un = Y^nMAM^n ' M aBAbY' (61) 

This means that for a local excitation, say on A {n = ua)-, the left eigenstate will take on 
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the form 

= J2^uK\(^o\ + E ^ujKiiKi (62) 

Ia IaIb 

where Yj^n = ^/^n- exact (full) bCC result this transforms into the product of the 

excited fragment-A state and dual ground-state of fragment B (see Sec. IV. E). 

B. Transition moments: truncation errors 

Now we are in the position to examine the truncation errors in the transition moments. 



Let us first consider the left transition moments (Eq. 45), 

r« = (vl/W|D|vl/-) (63) 

which may be written more explicitly as scalar products 

T(0 = yt ^(0 (64) 

of the left eigenvector y„ and a vector F*-'-* of basis set transition moments, 

= {^i\D\%') (65) 

associated with the biorthogonal (left) basis states. The basis set transition moments, being 
part of the general bCC representation of the one-particle operator D, exhibit canonical 
order relations 

0[Z«]=/.-l (66) 
as shown in Fig. 5b. A proof of these order relations is given in App. A. The scalar product 



(64) combines the left eigenvector components with the left basis set transition moments. 
This means that the (lowest) PT order associated with a specific class fi of eigenvector 
components is given by 0[FJ,F^'^] = 0[F^] + 0[F^^^]. Accordingly, the truncation after 
class fx leads to an error of the order 0[Y^^^] + 0[fII^]. In the case of singly excited 
states, the Cl-type order relations of the left eigenvectors (Fig. 3b) lead to the TEO formula 



(53). However, so far we have disregarded the effect of truncations in the CC ground state. 
For example, there are first-order contributions in F^p (that is the p-/i-components of F*-'-*) 
associated with the 2p-2h cluster operator T2 in the CC expansion of \'^q'^)- In conjunction 
with y^, being of zeroth-order, this gives rise to a first-order truncation error in the left 
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transition moments if the ground-state CC expansion does not comprise the 2p-2h cluster 
operator. This can readily be generalized, noting that the cluster operators are of the 
PT order /i — 1. This means that in the class fi components of the F vector, F_^^\ there 
are /i-th order contributions arising from the T operators of class /i + 1. Depending on 
the respective lowest order in the eigenvector components of class fi, F^, this leads to an 
additional truncation error, the order of which is by 1 smaller than that arising from the 
eigenvector truncation. The overall truncation errors for singly excited states are given by 
the following formula: 

{l/i, a even 

(67) 
|/i — |, /i odd 

In a similar way, one may analyze the bCC truncation errors for doubly and higher excited 
states; general TEO formulas are given in App. B.2. 



Likewise, the right transition moments (Eq. 46) 



Xn{<ifo\D\%') + Y,FPXin (68) 



can be written as scalar products of the extended right eigenvectors X'^ (see Eq. 38 ) and a 
vector F*-*"^ of right basis set transition moments. 



fP = (^ol^l^?) (69) 



associated here with the CC states of right expansion manifold. In Eq. (68), the ground-state 
contribution (/ = 0) is written separately. 



Let us consider the second term on the rhs of Eq. (68), which will be seen to determine 



the overall truncation errors. Using the expansion (35) of the dual ground-state, this term 
can be written as 

J2 F'f^Xi^. = 5^ ($0 1 ^ I + t^DX^ (70) 

I I 

where D is the bCC representation of £), 

Dij = ($/|^|^!}) (71) 
The order relations of D are shown in Fig. 5 (see App. A for a proof). The summation 



in the first term on the rhs of Eq. (70) is restricted to p-h components, [/] = 1, because 
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($o|-D|\l//) = for [I] > 1, and, thus, does not cause a truncation error whatsoever. The 
second term can be written in the more exphcit form 



reflecting the underlying partitioning with respect to excitation classes. To determine the 
truncation errors (at the level /i) one has to analyze the PT orders of the contributions with 
K = /i + l,A</i + l and X = fi + l,n<fi + l. This is described in App. B, where general 
TEO formulas are derived. For singly excited states, [n] = 1, the TEOs of Yj^DX^^ given 



by Eq. (67) 



It remains to inspect the first term on the rhs of Eq. (68), being the product of the 
ground-state admixture coefficient, x„ = — l^jx^, and the ground-state expectation value 



{'$o\D\'$q'^). The truncation errors in the Xn coefficient have been specified in Eq. (B.37) of 
App. B.2, which become 

^ " /i even 



Ote{ij) 



j/i + ^, /i odd 



(72) 



in the case of singly excited states, [n] = 1. To determine the truncation errors in the factor 



q\D\^'{^) , we expand the dual ground state according to Eq. (35), which yields 



(73) 



Here the first term is of zeroth order, and does not induce a truncation error. The second 
term can be analyzed in a similar way as the left TMs above. Here the TEOs are given by 



H even 



'/i + ^, /i odd 



(74) 



Since for singly excited states {[n] = 1) x„ is (at least) of 2nd order, the overall truncation 



errors in the first term on the rhs of Eq. (68) are given by that of x„ (Eq. 72), exceeding 



those of the second term (Eq. 67) by one. We note that for singly excited states the TEOs 
in the left and right transition moments are the same (see Table 1). 
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C. Transition moments: separability 



The separability of the left transition moments (Eqs. 63 65) is easily established. Ac- 
cording to 

= {K\(^o\Da + Ds K> lO = (KIDaK) = Ffl^ (75) 

it is seen that for local excitations I a the basis set transition moments for fragment A are 
identical to the corresponding moments of the composite system. Moreover, the basis set 
transition moments for mixed excitations vanish, 

Ffli, = (76) 

which means that the non-separable part of the left eigenvector does not come into play at 
all. Accordingly, for a local (fragment-A) excitation, n = ra^, we may write 

^i'^ = E^/.n^r = ^r (77) 
Ia 

where Tn^^ is the left transition moment for fragment A. In deriving this result, we have 
used that the local components for the composite system eigenvectors are identical to the 
components of the corresponding fragment eigenvectors, Yj^n = ^/^n- 

Whereas the left transition moments are separable notwithstanding the non-separable left 
eigenvectors, the right transition moments (Eqs. 68p9 ), involving separable eigenvectors. 



prove not to be separable. The problem here arises from the use of the dual ground state 
(\l'o|, more precisely, from the fact that a factorization of the dual ground state, 

(Wo| = (<|(^f| (78) 

is attained only in the exact (full bCC) treatment. To better understand the problem let us 
first assume factorization of (\l/o| and inspect Tn^ for a local excitation, n = ua- 

T(^) = xMD\K) + E Kn (79) 

Ia 



(r) 

Here the separability properties (60) of the right eigenvector 2Ln have been used; Fj denote 



right) basis set transition moments (Eq. 69). The ground state component Xn of the right 



eigenvector in the first term on the rhs is always separable, since 

Xn = -YlX^ = -Y^^X^ = x^ (80) 
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By contrast, the ground-state expectation value 



(81) 



is a non-local quantity, involving both fragment A and B. It should be noted that the 
separation of the bCC ground-state expectation value into the sum of fragment expectation 
values not only holds for the exact (factorizing) dual ground state, but also for truncated 



expansions according to Eq. (84) below. (In this sense the bCC ground-state expectation 
values themselves are separable quantities.) Being a product of a local and a non-local 



factor, however, the first term on the rhs of Eq. (79) is not separable. This means that the 



-B, 



non-separable contribution x„(\E'q in the first term of Eq. (79) must be cancelled 

by a corresponding contribution in the second term, to be identified in the following. For a 



local configuration, / = Ia, the right basis set transition moments (Eq. 69) become 



(82) 



Obviously, Fj^J is not separable. While the first term on the rhs is the fragment-A transition 



moment, F^"^^ = (\I'q'|D^|\1/;^), the second term is a non-separable contribution involving 

(r) 

fragment B. In the full bCC treatment the two non-separable contributions in T„ cancel 
each other. This is readily seen by inserting Eq. (82) in Eq. (79) and using that {^q I'^f^) = 



Y^o and 



(83) 



Ia 

For truncated bCC expansions, on the other hand, the non-separable contributions will not 
compensate each other, giving rise to size-consistency errors in the computational results. 

The non-separability of the right transition moments can be further elaborated by in- 
specting the general form of the dual ground-state. 



J A Jb 



B I 
Jb\ 



yjajb 

JaJb 



^4) 



applying both to truncated and full expansions. Here the local expansion coefficients are 
separable, that is, Yj^o 



^jIo and Yj^q 



YjgQ- This follows from Eq. (|34|) and the 



separability properties of M and y_ (see below). Using the expansion (84), the right basis 
set transition moments take on the form 



Jb 



-B 



(85) 
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With the help of Eqs. (80) and (81) we finally arrive at the expression 



IaJb 



-B 



Xi.nY;^j,,o{^jjDs\^^) (86) 



for the right transition moments. Here Tn"^^ denotes the right transition moment for frag- 
ment A. The non-separable contributions are identified as the second and third term on 



the rhs of Eq. (86). These terms cancel each other if the non-local eigenvector components 



factorize,that is, 



Y, 



IaJb,0 



(87) 



This can be seen by recalling that X]^/a0^^a« = ^^^^ (\I^o I = ('^^1 + Z] ^jso('^Js I- 



should be clear, however, that the factorization (87) of the non-local eigenvector components 



is equivalent to the factorization (78) of the dual ground state (84), applying only to the 
full bCC expansion. 

It may be of interest to see how the factorization of the exact dual ground-state eigenvector 



components derives from the explicit expression (Eq. 34) 



1 n 



tT\/r-l 



As is readily established, the local contributions to y_ are separable, 
and the mixed components vanish, 



The inverse of the bCC secular matrix is given by 



(88) 



(89) 



^ - Pa,ab 

~ -^BB Qa,AB 



- ^AB,AB J 



(90) 



where 



Pa,ab = -M^\Ma,abM/s^^j, 
Qa,ab = -^bb^b,abM^^^j^b 



so that, according to Eq. (34), the non-local eigenvector components can be written as 



Yaba = {v\M-a\Ma,ab + v'bM-b\Mb,ab)M-j,'^^aj, 
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(91) 
(92) 

(93) 



To proceed the secular matrix blocks Mab,ab, Ma,ab, and Mb,ab have to be further 
evaluated. Using a somewhat symbolic but largely self-explanatory notation these blocks 
may be written as 



Mab,ab = 1b X Maa + 1a X Mbb (94) 
Ma,ab = 1a X Mb,ab = 1b X (95) 

Proceeding at this symbolic level the desired result is readily obtained as follows: 

Yab,o = V^b{M-^\ + MI^bWaI^ab 

= v!bMa\M^I{1b X Maa + 1a x Mbb)M^]^,^s 

= v\Ma\ X v'bMb^b 

= Yao X Zbo (96) 

In a more stringent manner, the preceding computation can be performed on the matrix- 
element level, that is, by explicitly expanding all matrix multiplications. Here, the symbolic 



treatment according to Eqs. (93 96) may serve as a guidance. 

Let us recall once again that a factorization of the dual ground state cannot be expected if 
the configuration space is truncated. For example, assume a configuration space extending 
through double excitations and let Ia and Jb denote double excitations on fragment A 



and B, respectively. Then the factorization according to Eq. (87) would require that the 
configuration space of the system as a whole comprises quadruple excitations of the type 
IaJb, which, however, are not available in the truncated configuration manifold. 



D. CCLR form of right transition moments 

The derivation of the excited state CC equations in the framework of the linear response 
theory leads to the following separable, if more involved expression for the right transition 
moment [7j: 

tW = (^o|[^, Cr^Wo) - J2(^o\[[H, Cr^miiM + u;^)jj{^j\D\%^) (97) 

Here 

Cn = XktiCk (98) 
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denotes an excitation operator associated with the the n-th (right) excited state: l^'i'^'') = 
Xnl^'o) + C'nl^o'^)- Since the CCLR derivation starts out from a separable expression for 
a time-dependent ground-state expectation value, one may expect that the separability 
properties will be maintained in the further development. Nevertheless, it is reassuring to 
see directly that the CCLR form of the right transition moments is separable [I9] . Moreover, 



one will expect that the ordinary bCC (68) and the CCLR (97) expressions, while being of 
quite different form, must somehow become equivalent in the exact (full) bCC treatment. 
The absolutely non-trivial proof of this equivalence has been accomplished by Koch et al. |19] . 
In the following we will briefly review the separability of the CCLR right transition moments. 
The equivalence of the two transition moment expressions is addressed in App. C. 

To show the separability of the right CCLR transition moments we will suppose the 



general expansion (84) of the dual ground state, which holds both for approximate and 
exact (full) bCC treatments. Let n be a local excitation on fragment A [n = ua)- According 
to the separability of the right eigenvector, Cn consists only of local excitation operators, 
Cn = Yl -^Ia^Ia- Now it is easy to see that the first (commutator) term on the rhs of 
Eq. (97) is separable. Since [D, Cj^] = [Da, Cj^], the commutator becomes a local (fragment- 
A) operator, say Oa- As a consequence, in the matrix element (\1/o|Oa|^o^) the fragment 
B and non-local (AB) contributions in the expansion of (\l/o| are projected out, that is, 
(\l'o|Oyi|\l/Q'^) = (\E'o I^aI^o )• Now let us consider the second term on the rhs of Eq. (97), 
involving a double summation running over generic configuration indices J, J. The double 
commutator, involving the fragment- A excitation operator Cn and excitation operators Cj, 
leads to a restriction upon the indices J: there are no (non-vanishing) contributions for 
I = Ib (fragment B excitations), since [H,Cig] = [HB,Cig] and [[i^s, C*/^], C*„^] = 0. But 
what about contributions associated with non-local configurations / = Iab, not excluded by 
the double commutator term? To proceed let us inspect the matrix elements of {M + uJn)~^ ■ 



According to the separability structure of this matrix (see Eq. 90) the only non- vanishing 
matrix elements of the type {Iab, J) are those where the J index is non-local too, J = Jab- 
However, the non-local left basis set transition moments, Fj^ = ($j|Z)|\I^q'^), appearing 



as factors on the rhs of Eq. (97), vanish for non-local (mixed) configurations Jab (see 



Eq. 76), which means that non-local configurations can be excluded in both the J and the 
/ summations. To conclude: for a local excitation n = ha, the double summation on the 



rhs of Eq. (97) runs only over local fragment A configurations I a, J a- With this restriction 
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it is readily established that the three ingredients on the rhs of Eq. (97), i.e., the double 
commutator matrix elements, the matrix inverse, and the left basis set transition moments, 
are separable: they give the same results irrespective of being computed for the entire system 
or for fragment A only. 

E. Excited state properties and transition moments 

So far we have discussed ground-to-excited state transition moments required to com- 
pute spectral intensities. Now we will turn to excited-state expectation values (properties) 
for physical quantities of interest, e.g., excited-state dipole moments, and, more generally, 
transition moments associated with transitions between two excited states. 

In bCC form the general expression, 

Tnm = (^nl^l^m) (99) 

for excited-state transition moments becomes 

T_ = (vl/«|D|vl/M) 

= Xm{^^^\D\%'') + YiDX^ (100) 



Here D denotes the bCC representation (71) of a given operator D. The order structure of 
D is shown in Fig. 5, supposing here that is a one-particle operator; for a proof of these 
order relations see App. A. 

The truncation errors of the Tnm matrix elements are governed by the 2nd term on the 



rhs of Eq. (100). The secondary role of the first term can be seen in a similar way as in the 
right transition moments discussed in Sec. IV. B. We here skip the corresponding analysis of 
the first term, noting only that the order relations of two constituents, that is, the ground- 
state admixture coefficient Xm and the left (ground-to-excited state) transition moment Tn\ 
have already been established in Sec. 4.B. The truncation errors associated with the second 
term, being of the form of a vector x matrix x vector product, can be derived from the order 
relations of D and the respective left and right eigenvectors, as described more detailed in 



App. B.2. In the case of singly excited states {[n] = [m] = 1), the general formula (B.36) in 
App. B.2 simplifies to 

i/i — 1, a even 

OTEi^i) = { ' (101) 

— 1, /i odd 
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As the comparison with Eq. (53) shows, the excited-state transition moments for singly 
excited states (and a one-particle transition operator) have larger truncation errors (i.e. 
lower TEOs) than the excitation energies and ground-to-excited state transition moments. 

To discuss the separability we consider once more local excitations on fragment A, that 
is, n = Ha, m = rriA- In the exact case, where the left and right excited states can be written 
as fragment state products, 

the excited-state transition moments take on the manifestly separable form, 

Tnm = T± + S^^e^^lDBl"^^) (102) 

where 

T± = {^^^^IDaI^^^'') (103) 

is the transition moment for fragment A. Note that for diagonal (property) matrix elements 
{n = m) there is a contribution {^q \Db\^o)) corresponding to the ground-state expectation 
value of D for fragment B. 



Now let us analyze Eq. (100) in the case of a truncated bCC representation. The first 



term on the rhs of Eq. (100) is separable, 

^mTji — (104) 
as has already been shown in the preceding subsection. In the second term. 

Km = XiDX^ (105) 

the separability structure of D comes into play. As in the case of the bCC secular matrix 
(see Section IV. A), the separability structure of D can easily be derived. The result is shown 
in Fig. 6. The relevant subblocks Daa and Dab,a are given by 

Daa = D% + 1a{^^\Db\'^o) (106) 
Dab,a = 1a X F^S (107) 
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Here = F}^'" is the vector of left basis set transition moments for fragment B, as in 



Eq. (75). Using these expressions, Eq. (105) becomes 



r 

nm 



Y^An^AAKAm + Y^AB AB ,A^Am 



where 



Y 



t 

ABm 



,AB) 



(108) 



(109) 



is the non-local part of the left eigenvector Y}^, as specified by Eq. (61). Let us now first 
consider the non-diagonal case, n ^ m, where the left and right eigenvectors are orthogonal. 
Using the result 

Y^^DAAK^m = Yi^D^Xi (110) 



for the first term on the rhs of Eq. (108), as well as Eq. (104), leads to the following 
expression: 



T 



+Y 



AB,n^AB,Ad^Am 



(Ill) 



This means that Tnm is not separable due to the second term on the rhs arising from the 
non-local left eigenvector contributions. In the exact (full bCC) treatment these non-local 
eigenvector components factorize according to 



(112) 



that is, they form a products of excited state and ground-state eigenvector components 
for fragment A and B, respectively. Then the non-separable term vanishes due to the 



orthogonality of the fragment eigenvectors, H^^X^ = and the form of Dab, a (Eq. 107). 
For the diagonal case, n = m, the result is 



(113) 



IaIb 



where again Eq. ( |107 ) has been used. Only upon factorization of the non-local eigenvector 



components (Eq. 112) the correct result of Eq. (102) is obtained 



Again we may perform a brief symbolic calculation to demonstrate the factorization (112) 



of the exact excited-state eigenvector components The starting point is Eq. (61), where we 
may replace Ma,ab according to 



Ma. 



AB 



1a X Vb 



-1a X Y^Bo^BB 



(114) 
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to give 



-An 



Here the general relation (34), specialized to fragment 5, 



(115) 



(116) 



has been used in Eq. (114). Since is an eigenvector of Maa it follows that 



Y\^{uJn - Mab,ab) = Y^aS^aa - Mab,ab) 



-Y\n X MbB 



(117) 
(118) 



and, as a consequence. 



Y 



AnV^n 



M 



AB,AB 



BB 



(119) 



Using the latter result in Eq. (115) gives 



Y' — Y' ^ Y' 

^AB,n — ^An ^ ^BO 



(120) 



The factorization of the non-local eigenvector components is of course equivalent to the 



factorization {¥n^\ = {^ll'^'K^l of the expansion (^62|. 



Jl)A 



B, 



As we have seen, the ordinary bCC expression (100) for the excited state-transition 



moments and properties is non-separable, which here is due to the non-separable components 
in the left excited-state eigenvectors. Again, the CCLR approach results in an alternative 
separable expression [7j, reading 



I,J 

(121) 

where Umn = — ^n- The separability of this form can be shown in the same manner as 
in the case of the ground-to-excited state transition moments (Sec. IV. D). The equivalence 



of the CCLR form (121) and the ordinary bCC expression of the excited-state transition 



moments is briefly addressed in App. C. 
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V. HERMITIAN INTERMEDIATE STATE REPRESENTATION 



The bCC representation is a mixed or hybrid representation made up from the CC states 
and the associated biorthogonal states. Whereas the (correlated excited) CC states are 
genuine intermediate states, being based on the exact ground state, the biorthogonal states 
are essentially of Cl-type, that is, excited HF configurations. As was analyzed in the pre- 
vious section, the use of the biorthogonal Cl-type states as the left expansion manifold 
downgrades to a certain extent the truncation errors and separability properties of the 
bCC computational schemes. For the purpose of comparison, we will briefly inspect the 
properties of a hermitian intermediate state representation (ISR), specifically the ADC-ISR 
approach [181 El ESI |26l 132] , in the following. 

As the bCC representation, the ADC-ISR approach starts from the correlated excited 



states (Eq. 15) 



\^^j)=Cj\^o) (122) 
where |\E'o) now refers to the normalized ground state rather than to the CC parametrization. 
The correlated excited states (CES) can then be transformed into orthonormal intermediate 
states (IS), 

m ^ 

via a (formal) Gram-Schmidt orthogonalization procedure, in which successively higher CES 
classes fi are orthogonalized with respect to the already constructed lower IS classes u = 
1, 2, . . . , yU — 1. Within a given excitation class, symmetric orthonormalization is adopted. 
All the states are explicitly orthogonalized to the exact ground state, forming a zeroth 
excitation class (yU = 0). As a result of this (so far purely formal) procedure, one obtains an 
orthonormal set of intermediate states {"^j), 

= 6jj (123) 

being, moreover, orthogonal to the exact ground-state, (\E'j|\E'o) = 0. 

Representing the (shifted) Hamiltonian H — Eq in terms of these intermediate states gives 
rise to a hermitian secular matrix M, 

Mjj = {^>i\H-Eo\^>j) (124) 

and the associated hermitian eigenvalue problem, 

MX = Xri, X"fX = 1 (125) 
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Here f2 is the diagonal matrix of excitation energies, ujn = En — Eq, and X denotes the 
matrix of (column) eigenvectors. The n-th excited state can be expanded as 

|^„) = 5^Xj„|^j) (126) 
J 

in terms of the intermediate states and the eigenvector components The transition 

moments take on the form 

J 

where 

Fj = {^j\D\^o) (128) 

are denoted as IS transition moments. 

To obtain practical computational schemes, the Gram-Schmidt procedure is used together 
with Rayleigh-Schrodinger (RS) perturbation theory for |\l/o) and Eq, generating explicit 
perturbation expansions for the secular matrix 

M = M(°) + M^^) + M(2) ^ _ _ _ ^^29) 
and the IS transition moments 

F = + + + . . . (130) 

By truncating the IS manifolds and the perturbation expansions for the secular matrix 
elements and IS transition moments in a systematic and consistent manner, one arrives at a 
hierarchy of ADC(r;,) approximations, where n indicates that both the energies and transition 
moments of the lowest excitation class (singly excited states) are treated consistently through 
order n. An alternative and, beyond 2nd order, preferable derivation of the ADC-ISR 
perturbation expansions is the original ADC formulation [211 [25] based on diagrammatic 
perturbation theory for the polarization propagator 

As a distinctive feature of the ADC-ISR, the canonical order structure [IB] applies to the 
entire secular matrix (see Fig. 7), 

Mjjr^om-m (131) 

and, as a consequence, also to the eigenvector matrix X: 

Xjj^O{\[I]-[J]\) (132) 
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In analogy to the last paragraph of App. B.l, one readily obtains the truncation error formula 

= 20[X^^,J = - [n] + 1), /i > [n] (133) 

for the excitation energies. Here [n] denotes the class of the excited state n (as established 
by the PT parentage), and /i specifies the truncation level of the ISR expansion manifold. 
In a similar way, the canonical order relations 

Fj ~ 0{[J] - 1) (134) 



for the IS transition moments lead (via Eq. 127) to the following expression for the truncation 
error in the transition moments: 

= + 0[X^^,J = 2/i - [n] + 1, /X > [n] (135) 

For the lowest excitation class of the singly excited states {[n] = 1), the truncation error is 
2fi, both for the excitation energies and the transition moments. In Table 1, the errors for 
the six lowest truncation levels are compared to the corresponding CI and bCC values. 
The treatment of excited state properties and transition moments, 

T^m = (^nl^l^m) = X^DX^ (136) 

is based on the ISR of a general one-particle operator, 

Dij = {^i\D\^j) (137) 



Like the secular matrix (Eq. 129), the IS property matrix D is subject to a perturbation 
expansion, 

D = D^""^ + D^^^ + D^""^ + ... (138) 
Here the "shifted canonical" order relations 

Djj^0m-[J]\-1), |[/]-[J]|>l (139) 

apply, reflecting that a one-particle operator can couple HP (zeroth order) excitations of 
successive excitation classes. The product = -DX„ modifies the canonical order relations 
of an ADC-ISR eigenvector accordingly, that is. 



~ 0{[I] - [n] - 1), [/] - [n] > 1 
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(140) 



This leads readily to the expression 

= 0[X^.+i,J + 0[Z^^,J = 2(/i - [n]) + 1, /i > [n] (141) 

for the truncation error of excited state property matrix elements T„„. 

The ADC-ISR secular matrix is fully separable [T8| it has a diagonal partitioning 
structure as shown in Fig. 8, and the diagonal blocks are identical with the corresponding 
fragment secular matrices, that is, Maa = M^'^\ For a local excitation, say on fragment A, 
the fragment and entire system treatments give the same excitation energy, and the fragment 
eigenvector is part of the entire system eigenvector, in which all non-local components vanish 
{Xjg^n = Xj^B,n = 0). The IS transition moments are separable as well, Fj^ = fI^\ 
Together with the separable eigenvectors this ensures size-consistent results for the ADC- 
ISR transition moments. The case of the excited state properties and transition moments 
has been discussed in Ref. [32] . 

VI. CONCLUDING REMARKS 

The basic concept underlying the CCLR and EOM-CC methods for electronic excitation 
in atoms and molecules consists in a specific biorthogonal representation of the (shifted) 
Hamiltonian in terms of two distinct sets of states: on the one hand, the set of excited 
CC states based the CC ground state, and, on the other hand, the set of their biorthogonal 
counterparts. This results in a non-hermitian (bCC) secular matrix. The excitation energies 
are obtained as the eigenvalues of the bCC secular matrix, while both the left and right 
eigenvectors enter the calculation of the spectral intensities. The two sets of states are of 
quite different quality. While the latter (biorthogonal) basis states are essentially of Cl-type, 
that is, excited Hartree-Fock (HF) configurations, the excited CC states, by contrast, are 
formed by applying physical excitation operators to the exact A^-electron ground state. The 
resulting correlated excited (CE) states are expected to be superior to the simple CI states 
as they already account for a major part of electron correlation. The intuitive idea here is 
that electron correlation in excited states should not be completely different from that in 
the ground state. In fact, the use of CE intermediate states warrants distinctive advantages 
over the simple CI treatment. Foremost, this concerns the truncation error associated with 
limited expansion manifolds. To attain comparable accuracy, the manifold of CE states can 
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be truncated at distinctly smaller excitation levels than the CI expansions. Moreover, the 
excited CC states are intrinsically separable with regard to (hypothetical) non-interacting 
fragments, which, in contrast to CI, allows to devise size-consistent approximation schemes 
based on truncated expansion manifolds. 

However, due to the equitable use of the biorthogonal set of states being essentially of 
Cl-type, the bCC representation must be viewed as a CI-ISR hybrid rather than a full ISR 
approach. This becomes strikingly manifest in the split order structure of the bCC secular 
matrix, being canonical and of Cl-type in the LL and UR parts, respectively. As a result, 
the truncation error and separability properties are clearly superior to those of CI, but also 
weaker than those of a full ISR method such as the ADC-ISR presented in Sec. V. This is 
reflected in the general truncation error order (TEO) formulas derived here for excitation 
energies, transition moments, and property matrix elements. In the case of single excitations 
at truncation level 2 (that is, neglecting triply and higher excited configurations), the TEOs 
in the excitation energies and transition moments are 2, 3, and 4 for CI, bCC, and ADC, 
respectively. At higher truncation levels the gap between CI and bCC, as well as that 
between bCC and ADC widens. At the (already somewhat academic) truncation level 6 the 
respective TEOs are 6, 9, and 12. Of course, a given approximation may not exhaust the 
margin afforded by the respective TEO. For example, the error in the transition moments 
of the CC2 scheme is of PT order 2 (due to the first-order approximation used for the T2 
amphtudes), while the truncation error is of PT order 3. The ADC(2) approximation allows 
for a consistent treatment of (single excitation) energies and transition moments through 
2nd order, the TEO being 4. Let us note that the PT order of the overall error is only an 
indicator for the quality of an approximation scheme. A large error order does not by itself 
guarantee accurate results, but rather must be seen as a necessary condition for accuracy: 
a certain accuracy level can only be attained in compliance with a corresponding PT order 
of the characteristic error. 

The hybrid character of the bCC representation is also reflected in the separability prop- 
erties. The excitation energies, given as the roots of the characteristic polynomial, are sepa- 
rable, which, in principle, ensures size-consistent results at approximative levels beneath the 
full bCC treatment. It should be noted, however, that the separability of the eigenvalues 
applies strictly speaking only to the (fictitious) separated fragment model. Allowing for a 
small interaction between the fragments A and B, which is a more realistic simulation of an 
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extended system, the separability becomes blurred because the coupling block Mab,a (see 
Fig. 4) no longer vanishes, and local and non-local excitations may mix. What is problematic 
here is that the small matrix elements of Mab,a, associated with the weak physical cou- 
pling of the fragments, form products, e.g. in the characteristical polynomial, with the large 
(non-local coupling) matrix elements of Ma,ab- The mixing can become substantial when 
local and non-local excitations are nearly degenerate. This problem has been addressed by 
Helgaker et al. (see p. 684), but a thorough dedicated study seems to be still outstanding. 

The bCC eigenvectors do not perform uniformly, as the left and right eigenvectors are 
non-seperable and seperable, respectively. In the associated left and (ordinary) right tran- 
sition moments, both of which arc needed for the computation of spectral intensities, the 
separability properties are reversed: separable left TMs, as opposed to non-separable right 
TMs. This reversal is due to the use of the dual ground-state and the CC ground-state in the 
right and left TMs, respectively. As a result, the spectral intensities based on the ordinary 
bCC TMs are not size-consistent. This problem does not arise within the CCLR framework, 
as here a separable, though more involved expression for the right TM is used. In the full 
bCC limit the CCLR expression and the ordinary one, used in the EOM-CC methods, are 
equivalent. At approximate levels, however, the consistency of the left (ordinary) and right 
(CCLR) TMs may become an issue. At the simple CCS (singles) level, for example, both 
the left and ordinary right TMs are consistent through zeroth order (due to the use of the 
HF ground state). The right TM in the CCLR formulation, however, is consistent through 
first order. This means that the CCLR spectral intensity expression combines a zeroth-order 
left TM and a first-order right TM. One might be inclined to see this as an improvement, 
but it should be noted that a result comprising incomplete first-order terms may be inferior 
to that of a consistent zeroth-order approximation. The consistency problem emerges also 
in the finding that for transitions with low spectral strengths the signs of the left and right 
CCLR transition moments may differ, leading to unphysical (i.e. negative) intensities. This 
is to say that the CCLR results are not necessarily more accurate than those of EOM-CC 
as long as size-consistency does not play a role. A conclusive comparative test of EOM-CC 
and CCLR intensity results for smaller and medium-sized molecules would be highly desir- 
able. Such a study should also comprise excited state properties for which both ordinary 
non-separable bCC or separable CCLR expressions are available. 

Of course, the most obvious drawback of the bCC methods is the non-hermiticity of the 
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respective secular matrices M. As already discussed, both the right and left eigenvectors 
are needed when spectral intensities and property matrix elements are to be determined, 
which requires an additional effort compared to the hermitian eigenvalue problem of a full 
ISR method. In the case of degenerate eigenvalues, care has to be taken to ensure the proper 
biorthogonalization of the associated sets of left and right eigenvectors. More disturbing is 
the possibility of complex eigenvalues. Even though the underlying Hamiltonian is hermitian, 
so that the excitation energies obtained as eigenvalues of M must ultimately (in the full 
bCC limit) be real quantities, one may encounter problems in actual computations. As first 
noticed and analyzed by Hattig [35] , complex eigenvalues can occur in the vicinity of conical 
intersections of two excited state energy surfaces (or hyper-surfaces). Kohn and Tajti|36j 
have developed some ideas of how to deal with that situation, but so far there remain open 
questions. 

The present analysis of the bCC methods for (neutral) electron excitations in an A^- 
electron system can readily be extended to the case of generalized excitations, such as (A^- 
1)- or (A^+l)-electron excitations used in the treatment of ionization or electron attachment 
processes, respectively. Corresponding CC methods have been referred to as IP-EOM-CC 
and EA-EOM-CC [371 138| l39| HQ]. In the case of (single) ionization the neutral operators 
^ have to be replaced by the manifold 

{Cj} = {ck;cickCi,k<l;...} (142) 

of physical Ih-, 2h-lp-, 3h-2p-,. . . operators. As in the case of neutral excitations, the succes- 
sive (A^-l)-electron excitation classes of fih-{fi-l)p excitations are labeled by yU = 1, 2, . . . ; a 
corresponding classification, [n] = 1,2, ... , applies to the cationic energy eigenstates 
indicating the respective PT parentage. In the (A^-l)-electron case (as in the other general- 
ized excitations) the complication arising from an admixture of the A^-electron ground state 



in the neutral excitations (Eqs. 38 41 ) does not apply, which somewhat simplifies the ionic 
bCC equations. With a few obvious adjustments, the discussion and the findings for the 
neutral excitations can readily be transferred to the case of (A^-l)-electron (and the other 
generalized) bCC schemes. It should be noted that here generalized transition moments 
Tn = (\l/o|-D|^I'^^^) come into play, defined with respect to a suitable electron removal (or 
attachment) operator of the type D = ^ dpC^ . This means that in the bCC representation 



Dij = {^i\D\'^j) to be used in the analogs of Eqs. (70[71 ), the states ($/| on the left side 
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of the matrix elements are the A^-electron biorthogonal states (Eq. 16 based on the neutral 
operators |4]) . In the discussion of ionic state properties and transition moments according 
to Sec. IV. E, again a particle-number conserving operator D is to be considered. The cor- 
responding bCC representation, Djj = is a pure (A^-l)-electron representation, 



where both ($/| and are based on the operators (142). 

The analysis given here of the EOM-CC and CCLR methods from the perspective of the 
bCC representation has shown decisive advantages over the conventional CI treatment, but 
also distinctly weaker TEO and separability properties than those of a full ISR approach such 
as the ADC-ISR. It should be noted, however, that the latter approach is manifestly based 
on perturbation theory for the secular matrix elements and effective transition coefficients, 
behaving essentially like the Rayleigh-Schrodinger (RS) PT expansions of the ground state 
energies and CI expansion coefficients. This means that both the ADC methods and ground- 
state PT have the same condition of applicability, namely a sufficiently large energy gap 
between the occupied and virtual HP orbital energies. When the energy gap becomes too 
small, for example at bond breaking nuclear conformations, PT based methods are bound 
to fail. The bCC quantities (secular matrix elements and basis set transition moments), 
on the other hand, are based on the T-amplitudes of the CC ground state, which can be 
determined in a completely non-perturbative way. Yet this edge over methods involving 
PT must be relativized, as the CC approach breaks down as well in situations where the 
ground-state is no longer adequately described by a dominant single reference configuration 
(see Bartlett and Musial [H], Sec. VI. C, and references therein). The reason is that the 
usual single-reference CC ansatz is ill-suited to deal with the so-called static correlation. 
As a remedy for this deficiency, much effort has been devoted to developing multi-reference 
(MR) CC schemes [111 gSl iU iSl SHI ISl (for a more complete list of references and 
an introduction into the vast field of MRCC methods the reader is referred to Sec. IX in 
the recent review article by Bartlett and Musial [H]). However, the MRCC approach to 
ground and excited states is far more complex than the single reference bCC representation 
considered here, and it has to be seen whether really effective computational schemes will 
emerge. 
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APPENDIX A: ORDER RELATIONS OF BCC REPRESENTATIONS 

A general proof of the canonical order relations in the lower left (LL) triangle of the bCC 
secular matrix can be found in Ref. [TH]. A brief review of the derivation of these order 
relations is given in the following. 

Let us first consider the simpler case of a one-particle operator reading in second- 
quantized notation 

D = '^dpgclcq (A.l) 

where dpg = {(f)p\d\(f)q) denote the one-particle matrix elements associated with D. The bCC 
representation of D, 

Dij = {^i\D\^]) 

= {^j\e~^De'^\^j) (A.2) 

was encountered in the treatment of transition moments and excited state properties, as 



discussed in Sees. IV. B and E (see Eq. 71) 



The bCC representation matrix D has an order structure associated with the partitioning 
according to excitation classes, as shown in Fig. 5. In the upper right (UR) triangle one finds 
the familiar CI structure for a one-particle operator. This result follows along the lines of 
the first paragraph in Sec. IV. A. In the lower left (LL) triangle the canonical order relations 

Du = 0{[I]-[J]-l), [I]>[J] (A.3) 

apply, which is to be shown in the following. 



The operator in the bCC matrix element (A.2) has a finite Baker-Hausdorff (BH) expan- 
sion, 

e-^De^ = D + [f),f] + ^[[D,f],f] (A.4) 
terminating here already after the double commutator term because 

f = J2tiCi (A.5) 
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consists of physical excitation operators only, and D has at most two unphysical operators. 
Let us now write the T operator according to 

f = Y,% (A.6) 

in terms of individual class operators T^, /i = 1, 2, . . . . The T-amplitudes, being themselves 
subject of a well-defined (diagrammatic) perturbation theory, exhibit the order relations (see 
Hubbard 03) 

f^~0(/i-l), /i>l (A.7) 
This means, for example, that the PT expansions of the T2 amplitudes, 

— ^ ] t abij C abij (^-8) 

begin in first order. The Ti amplitudes (/i = 1), being of 2nd order, are an exception 
reflecting BriouUin's theorem. 



What are the consequences of the expansion (A.4) and the order relations (A.7)? Since the 



BH expansion (A.4) begins with D, there will be non- vanishing zeroth-order contributions 
to Djj for [/] = [J] and [/] = [J] + 1. Now suppose that / and J differ by more than one 
class, that is, [/] > [J] + 2. In that case non-vanishing contributions in Djj will arise only if 
there are terms in the BH expansion that are at least of rank r = [I] — [J]. Here, the rank 
of an operator is the number of its (or c) factors. For example, D is of rank 1 and the 
operators are of rank fi. Now it is readily established that the commutators [D, T^] and 
[[£), Tfj],T^] are of rank /i and fi + u—l, respectively. (A commutator of two operators A and 
B with definite ranks, a and b, respectively, is of rank a + b — 1.) To determine the lowest 
(non-vanishing) PT contribution to the Djj matrix elements, one has to inspect the terms 



of the BH expansion (A.4) having rank r = [I] — [J] (which is the lowest rank allowing for 
non-vanishing matrix elements) and find the lowest PT order of those terms. For example, 
[D,T2] is of rank 2 and PT order 1, which means that for [/] = [J] + 2 the PT order of 
Dfj is 1. In the general case, [/] = [J] + fi, fi > 3, terms with the required rank r = /i 
and lowest PT order are due to the [D,T^] commutators, being of rank fi and PT order 
/i — 1. Likewise, also the double commutator [[£), T2], T^_i] gives rise to terms with rank n 
and order yU — 1, but there are no rank fi terms with PT order lower than fi—1. This proves 



the order relations (A.3). 
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Let us note that the order relations ($/|£'|\1/q'^) ~ 0{[I] — 1) for the left basis state 



transition moments (Eq. 65) follow as special case ([J] = 0). 



In a similar way the canonical order relations 

Mjj = 0{[I]-[J]), [I]>[J] (A.9) 
for the bCC secular matrix (LL triangle) elements, 

= {^i\e-^[H,Cj]e^\<^o) (A.IO) 

can be established. Now we have to consider the BH expansion involving the commutator 
Kj = [H, Cj] and check the emerging transition matrix elements of the type ($/|0|$o)- Iii 
contrast to the case of the transition operator considered above, Kj is itself of PT order 1 
and of rank [J] + 1 (regarding here only the relevant two-particle part of the Hamiltonian) . 
The BH expansion 

e-^kj e^ = kj + [kj, f] + i [[kj, f],f] + l [[[kj, f ], f ], f ] (A.U) 

terminates after the triple commutator, since Kj has not more than three unphysical c^(c) 
operators. Let us consider a secular matrix element Mjj, where [/] = [J] + /i, /i > 1. 



Obviously, terms of the BH expansion (A. 11) do not contribute to M/j if their rank is 



smaller than [/]. As above, we may analyze the terms of rank r = [J] + fi with respect to 



their PT order. For yU = 1, the first term Kj on the rhs of Eq. (A. 11) is of rank [J] + 1 
and order 1, thus giving rise to a first order contribution to Mjj. For higher values of fi, it 
suffices to consider the commutators [kj,T^], being of the required rank r = [J] + fi and 
PT order fi. Again it is readily established that there are no rank r = [J]+ terms of lower 
PT order. 



APPENDIX B: ORDER RELATIONS OF CI AND BCC EIGENVECTOR MA- 
TRICES 

The order structure of the CI and bCC secular matrices give rise to specific order relations 
for the eigenvector matrices, which, in turn, imply the respective truncation errors in the 
excitation energies and transition moments. In the following we will first consider the CI 
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eigenvector matrix (B.l), and then turn to the left and right eigenvector matrices associated 
with the bCC representation (B.2). In the third subsection B.3 we shall show how order 
relations established only for a triangular part of a matrix can be extended to the entire 
matrix as a consequence of unitarity. 



1. CI eigenvector matrix 



The order relations of the CI eigenvector matrix rely on perturbation theory (PT) for the 
exact states. Let us first consider the familiar case of the ground-state, where the well-known 
Rayleigh-Schrodinger PT can be cast in the compact form 



|*o) = I'^'o) + 



(0) 



(B.l) 



Here, the usual M0ller-Plesset decomposition of the Hamiltonian, 



H = Ho + Hi 



(B.2) 



into an unperturbed (HF) part Hq and an interaction part Hj is supposed; |$o) is the (HF) 



ground-state of Hq with the energy Eq^\ and Qo = 1 — |$o)('^'o|- To determine the lowest 
(nonvanishing) PT order for a specific eigenvalue component. 



Xjo = ($j|*o) 



(B.3) 



one has to analyze the contributions arising from the expansion on the rhs of Eq. (B.l). In 



z/-th order the leading operator term is Hj. Due to the two-electron (Coulomb repulsion) 
part of Hi, the matrix element {^j\Hj\^q) vanishes if the excitation class of J exceeds the 
value 2z/. For the excitation classes [J] = 2u and [J] = 2u — 1, on the other hand, the 
matrix element gives rise to a non-vanishing z/-th order contribution. Obviously, there is no 
lower-order coupling between the HF ground-state and excitations of class 2i/ and 2z/ — 1. 
This means that Xjq is of PT order u for [J] = 2v and [J] = 2z/ — 1. 
This result can also be written in the form 



0[X 



11 even 
\{^ + 1), /i odd, > 1 



(B.4) 
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where // denotes collectively the configurations of class The p-h excitation class (// = 1) is 
an exception, as here Xjq ~ 0(2) due to BriouUin's theorem. The ground-state component 
Xqo is of course of zeroth order. In Fig. 1, the order structure of is depicted. 

Now we turn to the order relations of excited states l^n). Rather than using individual 
PT expansions, the following analysis will be based directly on the order structure of the 
CI secular matrix (Fig. 1). However, a remark concerning the significance of excited-state 
perturbation theory is appropriate. As is well-known, PT expansions for excited states and 
excited state energies are of little practical use because the possibility of small or vanishing 
denominators ( "dangerous denominators" ) in the PT expansions prevents meaningful com- 
putational results. In a formal sense, however, excited-state PT expansions can be generated 
analogously to the ground-state case, which then can be used to analyze, e.g., truncation 
errors of excited-state energies and transition moments. Underlying such a formal PT is the 
concept that each excited state is related to a specific CI state, 

from which it emerges when the scaled interaction, e.g., in the form \Hi, is gradually 
increased from A = to 1. For our purpose we do not need the individual PT descent 
of an exact excited state. It suffices to suppose that the exact states can be classified 
according their derivation from the (unperturbed) CI excitation classes, p-h, 2p-2h, . . . , etc. 
Analogously to the notation used for the CI excitation classes, we will denote by [n] the class 
of the exact state |\E'n), that is, [n] = if the excited state n derives from the fip-fih class of 
CI states. The classification of both the CI and the exact states allows one to partition the 
CI eigenvector matrix X into subblocks Xj^^, where jj, and z/ refer to the component and 
state classes, respectively. Fig. 9 shows the partitioning and the associated order structure 
of X. A general expression for the order structure is as follows: 

— a — u even 

II// — i/| |, iJ, — v odd 

The order relations of the eigenvector matrix refiect the underlying order structure of 
the CI secular matrix. We begin by considering the class of singly excited states {[n] — 1). 
Rather than dealing with individual eigenvectors, we can treat the entire set of class- 1 
eigenvectors at the same time. Let us therefore denote by the rectangular matrix formed 
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by all (column) eigenvectors X„ with [n] = 1. The eigenvalue equations for the eigenvectors 
of class 1 can be written compactly as 



HX^ = X^fli (B.6) 
where fli denotes the diagonal matrix of the p-h energy eigenvalues. Since any eigenvalue 



uin has an orbital energy (zeroth-order) contribution, that is, fii ~ 0(0), Eq. (B.6) leads to 
the following order equation 

0[HX,] = 0[X,] (B.7) 

This equation can be used to establish successively the individual orders of the component 
blocks, Xk,!- Obviously, the starting point is given by ~ 1 + 0(1), which merely reflects 
the fact that the singly excited states derive from the p-h CI configurations. To proceed 
we inspect the matrix-vector block products J2j H^jXj^i for successive values of the (row) 
index k. (To visualize these products it is helpful to write the order structure of H (Fig. 1) 
alongside the 2Li column matrix and fill in the successively determined order entries here, 
starting with the entry in the Xi^i subblock.) For k = 2 we may readily conclude that 

X2,i ~ Yl -^2,,^,-i ~ 0(1) (B.8) 

i>i 

where the first-order behaviour comes from the first term in the sum, i?2,i^i,i ~ H2,i ~ 
0(1). (Note that the diagonal eigenvector block behaves as = 1 + 0(1).) In a similar 
way, we may establish that also the next subblock is of first order, 2L31 ~ 0(1). For the 
4-th subblock the situation changes since the H^^i matrix block vanishes so that here and 
beyond the zeroth-order Xi i block drops out. One here obtains 

J]/f4,,X,-i~0(2) (B.9) 

where the second-order behaviour of X 4^1 derives from the first two summands involving the 
two first-order subblocks X2,i and X31. In the next step, one first-order block drops out of 
the matrix- vector product, but the second one, X^ i, combined with the -f^s^s block of the 
secular matrix, again leads to 2nd-order behavior of X^^i. Only for k = 6 the order jumps 
to 3, since now due to the structure of the secular matrix the first order blocks (and the 
zeroth-order block) no longer contribute to the matrix- vector product. Continuing in this 
way, the order relations of successively higher-class subblocks can be obtained. The general 
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pattern is that for each even k the order of the X^^i block increases by 1. Let us note that 
the procedure can readily be cast in the formally correct form of induction. 

Now we may proceed to the next higher class of 2p-2h states. Let X2 denote the 2p-2h 
eigenvector matrix with the (sub)blocks Xk,2- In accord with the PT origin of the 2p-2h 
states here the k = 2 (diagonal) block is of zcroth order, more specifically, X2.2 ~ 1 + 0(1). 
But the order of the first block, Xi^2 is fixed as well. As will be demonstrated in subsection 
B.3, the orthogonality of the p-h and 2p-2h eigenvectors requires that Xi^2 is of first order. 
In a similar way as above, the orders of the higher k blocks can now be derived successively 
from the matrix-vector products 

~ ^ 2 (B.IO) 

i>i 

It is readily established that the ^3^2 and ^4^2 blocks are of 1-st order, followed by two 
blocks of 2-nd order, and so forth. 

The procedure outlined here for the p-h and 2p-2h states can easily be extended to higher 
excitation classes, //. The diagonal block, X^^^ ~ 1 -|- 0(1) is always of zeroth order, while 
the order relations for blocks above the diagonal, k < fi, are determined by the orthogonality 
between the eigenvectors of class n and those of the lower classes (see subsection B.3 below). 
In the 3p-3h states, for example, the orthogonality constraint with respect to the p-h and 
2p-2h states requires the Xi^3 and ^2,3 blocks to be of 1st order. The full procedure 
for establishing the order relations of X, ascending both to higher excitation classes and 
higher blocks within a given excitation class, can, of course, be reformulated in a formally 
satisfactory way making use of induction. 

Once the order structure of the eigenvectors has been established, it is straightforward to 
analyze the truncation errors of the excitation energies. For this purpose one has to express 
an eigenvalue according to 

En = XiHX^ (B.ll) 

as an energy expectation value, involving the full secular matrix H and the exact eigenvector 
X_n- This expectation value can be written more explicitly as 

En = J2^LH.,xX,^ (B.12) 

k;,A 

where the greek subscripts are referring to excitation classes rather than to individual config- 
urations. To specify the error arising from truncating the CI manifold after class /i, one has 
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to inspect the PT order of the terms with k, = fi + 1, X < fi + l (or A = /i + l,K</i+l). Due 
to the structure of H it suffices to consider the diagonal contribution X}^_^_i^^H ^^i^^^iJ^^j^^^^. 
Since H 



0(0) the latter term is of the order 



(B.13) 



Using Eq. (B.5) this translates into the general truncation error formula 



fi — [n] + 2, ji — [n] even 



(B.14) 



/i — [n] + 1, n — [n] odd 

for the CI excitation energies of class [n], where, of course, /i > [n] is supposed. In a similar 
way one may analyze the expression (|9]) for the transition moments, yielding the following 
truncation error formula: 



fx 



n 



\n\ even 



/i — ^[n] + ^, [n] odd 



(B.15) 



Note that for the case of single excitations, [n] = 1, the general formulas (B.14, B.15) reduce 
to the simple expressions ^ and (11), respectively, given in Sec. 2. 
For completeness, let us note that the simple expression 



= /i - M + 1 

applies to the truncation error of property matrix elements. 



Tnn 



(B.16) 



(B.17) 



Here D is the CI representation ( 10 ) of a one-particle operator 



2. bCC eigenvector matrices 

The order structures of the right and left bCC eigenvector matrices X and Y, respectively, 
are shown in Fig. 10. (Note that X now denotes the right bCC eigenvector matrix rather 
than the CI eigenvector matrix as in the preceeding subsection B.l.) The bCC eigenvector 
matrices display both canonical and Cl-type order structures. The LL part of X and the 
UR part of Y are canonical, 

0[X^,]=0[Y,^]= fi-i^, fi>v (B.18) 
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while the UR part of X and the LL part of Y are of Cl-type, 

{Uu — u), a — V even, > 

(B.19) 
— u) + fi — u odd, > 

The order structures of the bCC eigenvector matrices can be deduced from the order struc- 
ture of the bCC secular matrix by an obvious generalization of the procedure used for the 
CI eigenvector matrix above (Sec. Bl). Here it is important that at each successive step 
associated with an excitation class /i both the right and left eigenvector matrices X ^ and 
Y must be treated in parallel, while the biorthogonality to the respective eigenvectors of 
the lower excitation classes, 1, . . . , /i — 1, must successively be taken into account. A formally 
correct way of translating the structure of the bCC secular matrix into the order structures 
of the eigenvector matrices can readily be spelled out. 

However, one may take also an alternative route, in which the order structures of the UR 
parts of both X and Y are determined directly. The order relations of the complementary 
LL parts are then obtained as a result of the biorthogonality, Y^X = 1. Let us first consider 
an eigenvector component 

= ($,|vI/M), [n]>[J] (B.20) 

from an UR block of X. Since {"^n"^) differs from the (normalized) exact eigenstate |\1'„) 
only by a normalization constant {Nn ~ 1 + 0(2)), we may consider ($j|\E'„) rather than 
Xjn, that is, 

Xjr, ~ ($j|^„,) (B.21) 



Using the expansion (30) for it is apparent that for [n] > [J] the bCC eigenvecor 



component is of the same order as the corresponding CI eigenvector component: 

0[($j|vl>„)] = 0[($j|x[/„)] (B.22) 
This establishes the Cl-type order relations for the UR part of X. 



In analogy to Eq. (B.21), the left bCC eigenvector components, 

F;„ = (vl>«|vl/0) (B.23) 
can be related to the normalized eigenstates, 

Yjn ~ (vl/„|^°) = i^nlCjin) ~ {^n\Cj\^o) (B.24) 
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where we have used the fact that the CC and the normahzed ground-state differ by a 
normahzation constant of the order 1 + 0(2). So it remains to show the order relations 

{^n\Cj\^o) ~ 0{[n] - [J]), [n] > [J] (B.25) 

for the "generahzed transition moments" (GTM) (\l/n|C*j|\t'o)- Let us first note that these 
GTM order relations are non-trivial. For a triply excited state {[n] = 3) and a p-h excitation 
operator ([J] = 1), for example, the (canonical) order is 2, rather than 1, as one might 
expect in view of two individual first-order contributions associated with |\E'^i^'*) and I^E'o^''), 
respectively. It is instructive to verify that these two first-order contributions, in fact, cancel 
each other. 

To prove the GTM order relations we may rely on the canonical order relations 

{^n\^j)^0{\[n]-[J]\) (B.26) 



as established for the ISR eigenvector components (see Eq. 132) discussed in Sec. V. Since 
the intermediate states {"^k) (including I^E'o)) form a complete and orthonormal set of states 
they can be used to expand the GTM (\['„|C'j|\E'o) as follows 

{^n\Cj\^0) = J2{^n\^K){^K\Cj\^o) 
K 

= {'^n\^K){^K\Cj\^o) (B.27) 

[K]<[J] 

where the second line is due to the fact that by construction the intermediate states {"^k) 
are orthogonal to all states CjI^o) with [J] < [K], that is, {^kICjI'^q) = for [J] < [K]. 
Now for [n] > [J], the order of {^nl'^x) decreases with increasing [K] and so does the order 
of kICjI'^o) ■ This means that the minimal order contributions in the sum on the rhs of 



Eq. (B.27) are due to those where [K] = [J]: 

(^„|C'j|^0) ~ 0[{^n\^j){^j\Cj\^o)] = 0[{^n\^j)] (B.28) 

The last equality follows from the observation that (\E'j|C*j|\E'o) is of zeroth order. This 
completes the proof of the canonical order relations for the eigenvector components in the 
UR blocks of Y. 

The biorthogonality of the left and right bCC eigenvector matrices exacts the order 
relations of the LL blocks of the X and Y matrices, as will be shown in the ensuing Section 
B.3. 
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Like in the last paragraph of Section B.l, the truncation error of the bCC excitation 
energies for general excitation classes [n] can be deduced from the eigenstate order relations. 
The starting point is the expression 



(B.29) 



where, analogously to Eq. (B.12), the second equation reflects the partitioning of the energy 



expectation value with respect to excitation classes. To determine the TEO at a given 
truncation level /i, one has to analyze the contributions where k = fi + l,[n] <A</i + l 
(set Si) and \ = fi + l,[n] < k < ji (set 82)- Here /i > [n] is supposed. The former set of 
contributions is given by 



Si = Y_ 



A=[n] 



(B.30) 



Due to the order relations in the LL parts of X (Eq. B.18) and M (Eq. 52), we find 



(B.31) 



irrespective of A. Together with the OR of (Eq. B.19) this leads to the following 

TEO formula 



;(/i - [n]) + 2, - [n] even 
'(/i - [n]) + |, /i - [n] odd 



(B.32) 



The 5*2 set consists only of two contributions. 



(B.33) 



because 7VfK,^+i = for k < /i — 1. The two involved subblocks of M are of first order, 
0[M"^^^+i] = 0[Af^_i^^+i] = 1, and the TEOs of the 5*2 contributions are seen to exceed 



those of 5*1. This means that Eq. (B.32) is the final expression for the truncation errors in 
the bCC excitation energies. 

In a similar way, one can derive general TEO formulas for the left and right transition 
moments, and the excited state properties. In case of the right transition moments, the dual 
ground-state eigenvector Fq comes into play. The Cl-type order relations of Y_q (see Fig. 3b) 



can readily be established by analyzing the eigenvalue equation Y}qM = —v^ (Eq. 34) in a 
similar way as in App. B.l, using here that only the p-h and 2p-2h components of are 
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non- vanishing, the latter being of first order. As to be expected, both the dual and the CI 



ground-state have the same order relations, as specified by Eq. (B.4) in App. B.l. 
The resulting TEO formulas are as follows. 

(i) left transition moments: 



l/i 



n 



fi — [n] even 
fi — [n] odd 



(B.34) 



(a) right transition moments: 



(Hi) property matrix elements: 



— [n] + 1, /i even 



- [n] + \, fi odd 



- [n]) + 1, fi- [n] even 

- [n]) + |, fi- [n] odd 



(B.35) 



(B.36) 



As above, [n] and /i denote the final state excitation class and the truncation level, respec- 
tively, where of course /i > [n]. 



The first expression (z) is obtained by generalizing the derivation of Eq. (67) in Sec. IV. B 



Here the error associated with the truncation of the excited state manifold follows from 



Eq. (64), using the order relations (66) of the left basis set transition moments together 



with the CI- type order relations (B.19) in the LL block of the left eigenvector matrix. 



In addition, one has to account for the error arising from the truncation of the ground 
state CC expansion, assuming here the same truncation levels in the ground and excited 
states. Supposing a sufficiently large or even complete ground state CC expansion, the TEO 



increases by 1 for even values of /i — [n] (first line on the rhs of Eq. B.34). 

In the case of the right transition moments (ii) the TEOs are determined by the second 



term on the rhs of Eq. (70). This term can be analyzed in an analogous way as the bCC 



eigenvalues by using the bCC order relations of the transition operator matrix D given in 
Fig. 5a rather than those of the bCC secular matrix. It should be noted that apart from 
the case of single excitations, [n] = 1, the right transition moments have larger truncation 
errors (lower orders) than the left ones. 
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In (m), finally, one has to analyze the the vector x matrix x vector-product in Eq. (100). 
The order relations coming here into play are those of D (Fig. 5a) and the LL parts of left 
and right eigenvector matrices. It should be noted that the expression [iii) applies also to 
inter-state transition moments for excited states of the same class, [n] = [m]. 

In the latter two cases, there is each an additional contribution, arising from the admixture 



of the ground state in the right eigenstate expansion (Eq. 41), specified by the respective 



(extended) eigenvector component x„ = —Y}qX_^ (Eq. 40). Using the order relations for Fq 
and the LL part of the right eigenvector matrix, one can readily derive the following TEO 
formula for the Xn components: 

OHW J^""'"'^'- (B.37) 

+ odd 

As discussed in Sec. 4.B and 4.E, the TEOs in the additional (ground-state admixture) 
terms exceed those of the respective main contribution. 



3. Order relations for biorthogonal matrices 

In the preceeding subsection B.2 the order relations for the bCC eigenvector matrices have 
been established only for the respective UR blocks, being of Cl-type in X and canonical in 
Y. Now we will show that the biorthogonality of X and Y requires canonical and CI- type 
behaviour in the LL blocks of X and Y, respectively. 

Let us first note that the biorthogonality relation 

Y^X = 1 (B.38) 

also implies that 

XY^ = 1 (B.39) 

which will be our starting point here. In this form the UR order relations of the first factor 
X (Cl-type) match the LL order relations of the second factor Y^^ (canonical), which for 
brevity will be denoted by Y' henceforth. For a graphical notion of the following procedure 
we recommend to place the X and Y' partitioning schemes next to each other and fill in 
successively the emerging order relation entries. 
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For the first row of X-blocks, Xi^k, k — 1,2,..., and the first column of V'-blocks, 
^M' k — 1,2, ... , the order relations are already given. In the second row of X-blocks and 
the second column of V'-blocks there is one undetermined block each, namely, X2,i and 
Y[ 2, respectively. The orthogonality of the second X row and the first Y' column can be 
expressed as follows: 

k 

What can be concluded from this with respect to the order of X2,i? Let focus on the first 
two terms in the sum, the remainder being (at least) of the order 2: 

X2,ir;,i + X2,2n,i + 0(2) = (B.41) 

Since the diagonal blocks of the eigenvector matrices behave as Y[ ^ = 1 + 0(1) and ^2,2 = 
1 + 0(1), respectively, and 1^2 ~ 0(1), the following relation holds through first order: 

X2,i + O(l)=0 (B.42) 

This means that X2,i must cancel a (non-vanishing) first-order contribution and, thus, is 
itself of the first order (more accurately: the lowest non- vanishing contribution in the PT 

expansion of ^2,1 is of first order). In a similar way, we may conclude Y\ 2 ~ 0(1), being 
a consequence of the orthogonality of the second column of l^'-blocks and the first row of 
X-blocks. 

After having completed the order relations in the second row and column of X and Y' , 
respectively, we may proceed to the third row of X-blocks. Here the orders of the first 
two blocks, X3 1 and 2, have to be derived, which in turn can be achieved by exploiting 

that this row is orthogonal to both the first and second column of l^'-blocks. Expanded 
explicitely through the first three terms these order relations read 

X3,in,l + ^3,2^,1 + ^3,3^,1 + 0(4) = (B.43) 
X3,ir;,2 + ^3,2^,2 + ^3,3^,2 + ^(3) - (B.44) 

The second equation allows us to determine the order of X3.2. Using Y'12 ~ 0(1), 1^32 ~ 
0(1), and 1^2,2 = 1 + 0(1), ^3^3 = 1 -|- 0(1) (as diagonal eigenvector blocks), we may 
conclude that the relation 

^3,2 + O(l)=0 (B.45) 
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holds through first order, which, as above, implies that ^3,2 is of first order. Using this 
result in the first orthogonality equation, together with Y[^^ = 1 + 0(1), ~ 0(1), and 
■^3,1 ~ ^(2)) yields through 2nd order 

X3,i + 0(2) = (B.46) 

which means that ^3,1 ~ 0(2), consistent with the canonical order relations. In a completely 
analogous way, the first two blocks in the third column of l^'-blocks can be treated, yielding 
the expected Cl-type order results, 1^'^ 3 ~ 0(1), 1^23 ~ 0(1) 

This brief demonstration has shown how the given Cl-type UR order relations of X and 
the canonical LL order relations of Y' (= Y^^) impose canonical order relations in the LL 
part of X and Cl-type order relations in the UR part Y^^ as a result of the biorthogonality 
of right and left bCC eigenvector matrices. Of course, this derivation can readily be cast 
into a formally correct proof by induction (see Ref. [IH])- 

In a related way, the unitarity of the CI eigenvector matrix (denoted X in App. B.l) 



can be used to extend the CI order relations (B.4), established in B.l only for the LL part 



to the entire matrix X. Writing the unitarity relation of the CI eigenvector matrix in the 
form X^X = 1, the order relations of the UR part of combine with those of the LL part 



of X = 1, similar to the product (B.39) of the bCC eigenvector matrices. The successive 
construction of the CI order relations in the UR part of X can be performed essentially as 
in the bCC case above. 

APPENDIX C: EQUIVALENCE OF CCLR AND ORDINARY BCC TRANSI- 
TION MOMENTS 

1. Right transition moments 



In the exact (full) bCC treatment the ordinary right bCC transition moment (Eq. 46 ) 



tM = i^olDCnlK) + x„(vI/o|^|vl>S'=) (C.l) 



and the separable CCLR expression (Eq. 97) 



TM = (^o|[Aa]|^S^)-5^(^o|[[^,0,],0„]|v]/-)(M + a^^^^^^^ (C.2) 
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are equivalent. Here, as in Sec. 4.D, 



(C.3) 



denotes the (right) excitation operator associated with the n-th excited state, \^n'^) = 
a^nl^o^) +C'n|^o'^)- This by no means obvious result was shown explicitly by Koch et al. [T9] . 
The following is essentially a reformulation and slight extension of the original proof, using 
the more transparent wave function notations promoted here. 



Let us start from the the ordinary bCC expression (C.l) and transform it successively 



into the CCLR expression (C.2). As a first step we make use of the commutator relation 



bCn = [D,Cn]+CnD, yielding 

= {^o\[D,Cn]\K) + {^o\CnD\K) + M^o\D\K) 



(C.4) 



To proceed let us consider the (trivial) identity 

{^0\CnD\K) = {^Q\Cn{H -Eo + UJn){H - + UJn)-'D\%') (C.5) 

and replace the inverse matrix operator on the rhs by its bCC representation. Noting that 
M' + Un is the bCC representation of H — Eq + Un, where M' is the extended bCC secular 



n —n 



matrix given by Eq. (27), the bCC representation of the inverse operator reads 

(M' + cu„,)-i 

where 



{M + UJnY^ 



(C.6) 



-u-WM + uJnV (C.7) 
This means that we can express {H — Eq + ujn)~^ as follows: 



Inserting this expansion into Eq. (C.5) yields 



i^olCnmo) = {^o\CniH -Eo + ^„)|^?) (M + ^„)7] {^j\mo) 

-Xn{{<l>0\D\%') +J2^nWnj{^j\m^)} (C.8) 
J 

In deriving this result we have used that {H — Eq + ojn) |^o^) = and (^o|C'n|^o'^) = 

—Xn- Note that at this point another a;„-term comes into play, augmenting the third term 
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on the rhs of Eq. (C.4). The reformulation of Eq. (C.5) is still not complete. To proceed 
the matrix elements (\[^o|C*„(i^^ — Eq + u;„)|\E'5) in the first term on the rhs of Eq. (C.8) can 



be expressed according to 



10 



(C.9) 



in terms of the double commutator [[H, Ci], C*„]- Here we have used the eigenvalue equation 
for the n-th excited state in the form 



(C.IO) 



Moreover, recall that (\I'o|C'/|^E'o'^) = ^/q and the operators Cn and Cj commute. Insert- 



ing Eq. (C.9) on the rhs of Eq. (C.8) constitutes the final step in the reformulation of 



{^Q\CnD\^!'^Q). Using this result in Eqs. (C.8) and (C.4), respectively, validates the second 
term in the CCLR expression ( |C.2[ ), but also introduces a third x„-term. 



(C.ll) 



I,J 



due to the second term on the rhs of Eq. (C.9). It remains to show that the three x„ terms 
cancel each other, that is, 



(C.12) 



where the contributions 1, 2+3, and 4 on the Ihs arise from Eqs. (C.4), (C.8), and (C.ll) 



respectively. The contributions 3 and 4 can be combined and further evaluated in a compact 
matrix notation as follows: 



uJn{w' + Yl){M + 



cj„(-cj- V(M + uj^y^ - v'M-\M + uJn)'^) (C.13) 

-ujnv\uj-^ + M-^){M + Lu^y' (C.14) 

-u;nv!u~^M-\M + Un){M + Un)-^ (C.15) 

-v'M-^ = yI (C.16) 



As a result the sum of the three Xn terms becomes 

where the cancellation now is obvious, as (\I'o| = ($o| + S^joi*^^!- 



(C.17) 
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2. Excited state transition moments 



In a similar way one may show the equivalence of the ordinary bCC and the CCLR 



expressions for excited state transition moments and properties, the latter reading (Eq. 121) 



= (<'^|[D,C7^]|vI/-)-^(vl/^'^|[[i7,C7,],(7J|vl/-^ 



7(0 



i,j 



where Umn = — ^n- As above, we start from the ordinary bCC expression ( |100[ ) 

T 

and use the commutator relation DCm = [D,Cm] + CmD yielding 



:vl/iVC'n.|^S'^)+X„(vl/iVl^S'^) 



Analogously to Eq. (C.5), we consider the identity 



7(0 



7(0 



(C.18) 



(C.19) 



(C.20) 



(C.21) 



and use the the bCC representation of [H — Eq + uj^n) ^ to further evaluate the rhs. Note 



that the only difference to Eqs. (C.6) and (C.7) is the replacement of Un by tOmn- This leads 
to 



7(0 



I,J 



i,j 



7(0, 



using here ICml^l'o^) = 5mn- Let us consider the 2nd term on the rhs, which is non- 
vanishing only for n = m due to the Kronecker symbol. Since Umn = for m = n, and 



(C.23) 



the latter term becomes Smn{'^o\D\'^Q'^), thus reproducing the 3rd term in the CCLR ex- 



pression (C.18). As in Eq. (C.9), we now may introduce a double commutator according 
to 

. . . . . . ^^^^^^ 



Here we have used the eigenvector equation for the m-th excited state in the form 



(C.25) 
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7(0 



and the relations |C/|\l/o^) = Yj^ for the left eigenvector components. Using this result 



in Eq. (C.22) gives 



+ u:^n)l] {^j\DM) (C.26) 



i,j 



The first term on the rhs is seen to reproduce, via Eq. (C.20), the 2nd term of the CCLR 



expression (C.18). It remains to inspect the 2nd term on the rhs, containing the factor Xr. 



Since y„ is a left eigenvector of M it follows that 



u ^Y^ 



(C.27) 



Thus 



X>^;n(^ + ^n.n)7](<f jI^I^S^) = (^iV I ^o'^) (C.28) 



which cancels the original Xm term on the rhs of Eq. (C.20). This concludes our proof. 
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TABLES 



TABLE L Truncation errors (PT order) for excitation energies, transition moments, and excited 
state properties of singly excited states: Comparison of CI, bCC, and ADC approaches at the 
lowest 6 truncation levels. 
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FIGURE CAPTIONS 

Fig. 1: (a) Order structure of the CI secular matrix H. The subblocks correspond to 
the partitioning with respect to excitation classes of iip-ixh (//-particle-//-hole) excitations, 
/i = 1,2,.... The entries 0, 1 indicate the "PT order" of the blocks, being here simply 
(in the diagonal blocks) or 1 (linear in the electron repulsion integrals); vanishing blocks are 
indicated by dashes, (b) Order relations of CI eigenvectors associated with singly excited 
states {Xph) and the ground state {Xq). The entries denote the (lowest) PT order of the 
respective eigenvector segments. 
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Fig. 2: (a) Block structure of the CI secular matrix H corresponding to the partitioning 
associated with the separate fragment model (see text) . (b) Structure of CI eigenvectors for 
a local excitation (on fragment A). 

Fig. 3: (a) Order relations of the bCC secular matrix M. As in Fig. 1, the block structure 
reflects the partitioning of M with respect to excitation classes n = 1,2, ... . The entries 
denote the (lowest) PT order of the matrix elements in the respective subblock; vanishing 
blocks are indicated by dashes, (b) Order relations of bCC eigenvectors: Xp^ and Yph 
denote right and left eigenvectors for singly excited states; Yq is the "dual" CC ground 
state. 

Fig. 4: (a) Block structure of the bCC secular matrix M corresponding to the partitioning 
associated with the separate fragment model (see text), (b) Structure of right and left bCC 
eigenvectors for a local excitation (on fragment A) . 

Fig. 5: (a) Order relations of the bCC representation £) of a one-particle operator D. Block 
structure and entries as in Figs. 1 and 3. (b) Order relations of the left and right basis set 
transition moments (see text). 

Fig. 6: (a) Block structure of the bCC representation £) of a general operator D with 
respect to the separate fragment model (as in Figs. 2 and 4). 

Fig. 7: (a) Order relations of the ADC-ISR secular matrix M. Block structure and entries 
as in Figs. 1 and 3. (b) Order relations of ADC-ISR eigenvectors for singly excited states. 
Fig. 8: (a) Block structure of the ADC-ISR secular matrix M with respect the separate 
fragment model, (b) Structure of ADC-ISR eigenvectors for a local excitation (on fragment 
A). 

Fig. 9: Order relations of the CI eigenvector matrix X. Block structure and entries as in 
Figs. 1 and 3. 

Fig. 10: Order relations of the right and left bCC eigenvector matrices X and Y, respectively. 
Block structure and entries as in Figs. 1 and 3. 
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Fig. 9 
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Fig. 10 
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